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• Rotational Equations of Motion 

• Structural Analysis - Boundary Conditions 

• Generic Modal Equations 
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• Consider SCOLE configuration without offset of the 

mast attachment to the reflector and without flexibility 

• Consider SCOLE configuration without mast flexibility 

but with offset in the direction of orbit (strawman) 

• Consider SCOLE configuration with offsets in two 

directions but neglecting mast flexibility 

• Consider general SCOLE system dynamics 

• IMPLICATIONS FOR CONTROL STRATEGIES 



Parametric Study of the System 

Let us assume that the interface point between the reflector and the 
mast is at the center of mass of the reflector 
♦ X - 0 ♦ A - 0 - C 5 - Cg 

Under this assumption, the equation becomes 

. <tyc. i /il Cy e , uaCOc*^) 

- C,/4C (I. -r„) = o (2.5) 

which in the absence of gravity gradient, yields the following first 
integral of the motion: 

- & +<Ve. a. Wc, **(at+4>) 

+ C4 4c, c«[*(n.-c + <j>)] = < (2 .6) 

This equation is plotted in the phase plane (9',e) for different values 
of u and Q. (Figs. 2.2) 

Floquet Analysis 

The angular motion about an axis perpendicular to the orbit plane in 
the absence of gravity gradient is described by: 

® = /JwiTI'C J & (2.7) 


2.3 
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Case 2. 


No gravity gradient, but offset. 


P(t) 


^ COT - tl 
1 C 1 


- — - ft sin At 
C 1 


CZ(.t)] - CP( T )] CZ(t)3 

^11 * p 11 z 11 ♦ p 12 z 21 O ) 

^12 m p 11 z 12 + p 12 z 22 ( 2 ) 

Z 21 " p 21 z 11 * p 22 z 21 (3) which becomes 

Z 21 " Z 11 since P 2 i - 1 and P 22 - 0 


z 22 ■ P 21 Z 12 + p 22 z 22 (**) which becomes 
z 22 ■ z 12 

from (3) Z 21 - i Uj substituted into (T) yields 
• • 

Z 21 * P 1 1 Z 21 + p 12 z 21 

similarly from (4) Z 22 - Z 12) substituted into (2) yields 

z 22 • P 1 1 z 22 * p 12 z 22 
C5 j ' 

since — - constant dt P 11 “ p 22 

Then 


*21 - P 1 1 Z 21 + P 1 1 Z 21 ’ Jt (P 1 1 z 21 ) 
and Z 22 - P n Z 22 - P n z 22 - ( p nZ 22 ) 

These two last equations are integrated and the following results 
for Z 2 i and Z 22 obtained. 

Z 21 " p i 1 Z 21 + Kt 
z 22 - P 1 1 z 22 + K 2 

but from (3), Z 21 (x) - Z n (T) and from (4) 
z 22< t > - Z U < T > 

Therefore, Z 21 (0) - Z^iO) - p n(0)Z 21 (0) + ^ -► - 1 

since Z n (0) - 1 and Z 21 (0) - 0 


Z 2 2<0) - Z 12 (0) - 0 - P n (0)Z 22 (0) ♦ K 2 

Cc Ca c * 

* " c* * — * K 2 * 0 or K 2 • - 

11 C 1 C 1 

The two last equations integrated once, yield 

z 21 - P n Z 21 ♦ 1 

• C 2" C 5 

z 22 “ Pl1 z 22 ~ H — “) 

C 1 


orsg’kal 

OF POOR QU3 


to 


^ Soluti on of the first order equations 

ir ~ p i i z 22 - - (~^-> 0) 

dz. 


e presence of g^p22 and PnZ 22 in the equation suggests a. product of the 
type *(t)Z 22 (t) 

but |p (*Z 22 ) ‘ J| Z 22 ♦ « |p Z 22 (2) 

Multiplying (1) by <Kt) yields 


♦ P 22 “ ♦ p 11 z 22 


which can become 


♦ <%?*) 

c i 


1 Q aQ 

JP (*Z 22 ) - - <P (-£~) 

C 1 


if one can find $(t) (integrating factor) such that 


d<fr 

dr “ ** * p 11 


in 


<Kt) - J-Pii dr - |-^2>cos flr dx + j ^ dr 

In <p ( r ) - - sin Qt 
c i n 


- ?x * K 

c X C 1 

or $ ( t ) - exp [--2- sin Qt] .e^f+K) 

c l ft C 1 

from d ( z 22<t ») - ( c2-cs ) 


dx 
( C5-C2 ) 
c. 


dx 


K] dz 


Z 22 ' f | • - 

Z 22 - exp [.Si. - Or r - K J ( C £zfi) £lZ + 

r _ T c, r \ c, 'J rL c *a . c, t 

•* c-s^ «i* •- sr r ^ - utf-^ i r 4* - - 

• Xp [*%] - 1 * C ‘%, * («» %,/•{ r - - • 

Therefore , exp fjak Qm XtZ + C* Jr U( Z + t .. - 

z 22 m exp 

^22^0) 


L^oC| ^ xCi i -* + “ 

1 ( c ^> ^( c +°sg*x! + +K) 

1 ♦ k, - 1 ♦ k, - 

C 1 1 1 c 5~ c 2 


2.11 


)? 


z 22 . «p[A»»Qt- |5tJ( i + + /c iz c t j i r% + 

a ‘"7 ■ z,2,t) ■ exp [& a tjncj.) , jW/c + .-...» ' 

<-f ^ - c_, cj (, : t jc^ c j r . .. j 

Z, 2 - exp [i(i« -«* C] ^,1.M t |(Vaot<,| ( cw.j 

* c^/jt *f( *JSS£^(' V*)' V 

Z 12 (t) - exp J**az~£!l-c]f£LbtU3az) + f&t-Ct){Cz tmJLV >Ct j ( 


Zi2(t) - exp 


f— faU nr » £s 

L-Q-c, c.) 


rJ[*^.(c»UT-oJ[i + ^s-c t r + ( c *~cij z Ci 


“d^ * p 11 z 21 * 1 where P n * ^ cos !Jt - ^ 

C 1 C]_ 

Integrating factor 4 ; . - ^»p 1t 

- expjl^^ar^^c^Kj and 

•«p [-£_ •* a« ♦ |r t] > i * (csj&*) c T ^ ~ • 

Integrating term by term yields, 

Z 21 “ exp — £iCj£‘C + (^rdZi) — + K. J J 

since Z 2 -|(0) - K’ - 0 

2 2 ,(t) - e* P [£i -g Z ][* + t C ^f*) % + j 

Zji • z,,(t) • exp | ,r ]{ ,+ “*• •• 1 
♦r* «p[$fe/* a * -%T][Z ♦ (C^.^ + f*z£<pl+ -J 

z n ( T ) - ex P [c^ Cl ^ar-c^,rj[i + (.Cs-c* t c*c»iii:-cr)r. + ... , 

Zii(t) - exp fifUA Q.Z — c% £j (C 03 ■flX-' |J "£ + Of ^ S^re-l) ; 

♦ J 

It is seen that Z n (0) - 1 


ORIG’ivAL PA^.2 «3 

OF POOR QUALITY 



a iUJb 








Fig. 2.2b SCOLE System-Phase 







FLOQUET STABILITY ANALYSIS 
2D SCOLE OPEN- LOOP SYSTEM 

Offset of the mast attachment point on the reflector 
results in an increase in the number of stable point 
for the lower frequencies 

Number of stable points increases for MR/Mm > 1.0 
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A. Angular Momentum of the Shuttle About ics Mass Center. G 


The angular momentum of the Shuttle, taken as a rigid body in a circular 
orbit, consists of contributions due to rotation about its center of mass 
plus the translation along the orbit. 


HS : 


where 


“ s/ Ro 


905,443 

0 


6,784,100 
145, 393 0 


in R(x,y,z) 


(1.9) 


145,393 

0 

7,086,601 


fry] 


( 1 . 10 ) 


us/ Ro " “* 1 + “y 1 + k 


( 1 . 11 ) 


B. Angular Momentum of the Beam about G 


Consider an element cf mass, dm, of the beam located at some point. 


p, such that GP - r + q - r 

o 


where: 


(1.12) r Q * -zk is the position vector of p in the undeformed state 

A A 

(1.13) q(z»t) - u i + V j in which, u and v are the x and y com- 
ponents of the mode shape vector. 


The angular momentum of dm about G, d Hm/ G is given by: 
d lte/ G - r x ^ (-Rk-+r) | RQ dm 

dita/ G - (?x %R i + rx ~?| Bo ) dm 


(1.14) 


< 


3.6 
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Which is expressed explicitly as: 

?fh **/)*, = A +(v+t$«+jc4,)f<, jt 

-) g/f** 

rx #/*o * 

«« r . ub stic utlng the differed teras lMO equaclm (1 . U)> the fo 
expression results: 

' — -•> 

= {[}( '*■ +j*ea*J£ 

}(ji -aJj v-)+ -a>*z>) _> Jj 

+[- MX , -U(V-*e^u). V(jI-uj 4 *u 4(juat a , s) 

Since u (2,t) - E P “ (t)s“(z) and v(z, t ) . E p n (t )s n (z), 
we consider for one mode in the open-loop situation. 


( 1 . 16 ) 


U - -u»; Sin («;t + a)s x (z) and 
"'S Sin (w y t + Y)s y (z) 


( 1 . 17 ) 


3.7 
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Assuming small elastic displacements such that, ^ X 


2 /_\ /„2 


and Sj. (*)/£ « 1, and dividing dff/ c by Q£ 2 , where Q is an assigned ' 

frequency and l a reference length, then. 


+ (- uj + Ctyjtt-J it jfe/} 

where p is the mass per unit length of the beam. After multiplying both 
sidesjDf this equation by Of 2 , there results: 

<=///*,/<; - + £j*&o -£U +£**%)£ 

+ (- vficDo + £M.Cd x +£2^0^) it J ( 1 . 19 ) 

The total angular momentum of the mast about G is obtained by inte- 
grating (1.19) over the total length of the mast, 

= J d/y*r/4 


(1.19) 


( 1 . 20 ) 


Th. e« terns appearing la dH./ 0 are Integrated using integral tables- e.g. 

r - S 4 *“(&■+?) J[ % * 4*^ 

-f *£+*(£**?} fa <s/&*)+ &/■&&. u&i / i 

A * ‘ + i? m J 

+ c *(- + 4 „ ^ue + j. )] 

* fit fi*> /if /9£/J 
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c * Angular Momentum of the Reflector about G. 

Since snail deflections for the beam are assumed, the reflector 

be assumed to be located at a constant distance from G, the Shuttle 
center. 


can 

mass 


Using the transfer theorem for the angular momentum, (See Appendix 


□XA) 


(i * 24) 


_ = Ir A cS r/Ko + HGZ £(6d,)j K 

where Ir/ 0l and I r / Ro - S t / s + »;/ 3 ( 1 . 25 ) are both expressed ln'the 
3a “ coordinate system, , moving wich the re£lector . In 

R 2 (principal axes of Inertia of the reflector) , 


Hr/ 


* Ir l 

0 

0 


4,969 


0 

0 

CM 

M 

0 

m 

0 

4,969 

0 

0 

» 

0 

Ir 3. 


„ 0 


9,938 


Mr/t, = ed.l +^+^4. + <l ! .t.i.0ty , + §r'6 l . 


( 1 . 2 ?) 


with j ’ » sin * i_+ C os $ 1 

r z r J 2 

therefore. 


^ r M> = + +9rA«t<$ r Z 1 + 8tB<$ r £i (1.28) 

$r ii 


(I. i6) 


3.10 
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( 1 . 34 ) 


« 


Hr/G - CD/ Jr, £ y. cD u Z z j\ + cojr, 


Where 



RfA-^yi-^v-c^yT, +XcO i Z,)-Kd>Y- eix)Zj- 


■jv- ^^^.^y^. tc4xZi , (aty ._ y 

' fa ~ **rZi * c4XZ t t (a,,y-al X ) 5} J 


3.13 


34 



D ‘ Aguiar Momentum of the System about G 


The angular momentum of the system - the sum of the angular momentum 
of each component evaluated at the same point 



— (oi Z+c$Zs)Z J/ -r (<$,&+ cqz,)£ 

+ f'flty <*>(*!*+'<)£- $^o4*+r)& - a*. £] 2 

3 

4.[cCoRP+ at! + a} } £]£ 

+[ak «)£ + u, cntyt+yjfi _ C&, Zen /$*+?) fa fej 

+ tyrf(tt+c(v+y))2 -(at+cCM+x))^ (b(*+x)-a(-<H-y ))£ / 

+ c45,Xt + aj 4*^ + a> s Zj £ (1M) 

In the expression for the totnl angular »o«eetu», the lest ter. will now 
be expressed In Bfx.y ,zi by slxply transferring i 2 , j' 2 , „ d ^ lnt0 ftmc . 
tions of i, j r and k as follows: 

•4 =• C*>&cn9r2 +M&cn</ f .f + / &u$ r /}^i l i, r £ 

+ (A m 4r&*> Sir &><£. f 



3.14 
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Rotational Equ ations of Motion. (Torque Free) 

In the linear range, 

Hx + /Yy - O (a) 

=> I ((/s~ + (0 + a) o pJf /3 + Mf[(<p'+O)o <t)cx> ((*£{+*]£ 

- C^i ( _ Ct)j*<un(tyt +jj£ - (<£- cOo fl)£ 3 J 

& 

+(£- cA 0+ £_ cAgJJ,, y. <u],g ZrJ+ (witp. a], 0) Z s , 

+(c*.<p + ag t)(r u . I 3t ) y - cA 0) cm(uiU«)£ 

V- (tit-la], i)e»(u$ti-]r)£ _ (ax>6 -w?)R antyt +r) fo] 

+ £&*%Z r -(%u£+eA0+ct£y.+ c*g r )Z s + Mr(u+X) fo-aljv 

+ J{cofj-a3.{)_ (uk0tu3g^)ji - X (tdfr+akg + (**?+&$) 
f coif r xJ- Mr Y) ha], 0- cot) R , (t>- w.)Z +( 

+ (lQ,<fi +a&pjif+ y(ci3o% + a).£ r + aioQ+alif’)- alZ&r Jcj 
+'Mr x ![#+( f- ah <p)e y X ($ + eap* £)] + MrXwX 

■+MrY[y~(ty- + ty.- CX>.(0-£) -}C(&+8r)J = O (IMS) 


3S 


3.18 



Hy + a £ H x - cO x - o (8j 

0 

-Ity (tf+rtopj - &£?+ cAfa-oO. yg Z 

+ 1(0- Br)- Mri[Il ~ 9i - V (0 +c<)op + $ r +a?<>%) 

+ COo &r Jcj -SMrXcOo u _ MrJCY( fi-COo <P+<fi r -u)o$ r ) 

- MrX*(B+&r) - M r y((/,-cOo<p)(eOoZ+a),t) = 0 ([-*}) 
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4 

■+ COk Hy - a)y/~/x - o (C) 

^ J( > i' + ((j>-LMo<yjZ ]3 4 - aie ft a )^ f 

+ C0z(y.-«3o<t>)4^(alt ia j^ i + fi cetyl y)£_ (ugfr-aJkjjtmtytty)^ 

^zy,^] + Hr(utx)[vt <yZ-a). <pt+ <i>4 
eo.9-M * X[<tKa.+j>+u.<f+g r -a>.<y r j] + 

x (?+'*’¥■*&)- MrjvtyjjA'-Se -<zSV- 
- y(ft 0)0 + $ r + &>¥+ <P) + ti>o X 9 r ] - Mr tffctitRi- u - (&-c0»X 

-($ + *.})*+ Uo&rX- yfr Uf +f\ ea ,p +$ r )7^.4) <A Zj t 

* ■p(9--‘*0)(‘**£ l + cdfrifa. &“■£) - cA(p_~e*.$)Z. i 

w 


4 

4 


- &* (9-t&4J+ Afy _ ea.^ % .(Ag fa , eofiZr 


+ Mrcdofve +(<£+0,0)4*+ Xi (£*.+$ 
+ 4+1*?)- cCo^Xt + YcOc u + y*fc+ifi r - to >j.ea.& ) 
-X(&-ctio+8r) — o (j-49) 


3.21 
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2. STRUCTURAL ANALYSIS 


A * Governing Differential Equations 

The governing partial differential equations for the system (beam) 
are comprised of two one-plane-bending equations (2.1) and (2.2) and one 
axial torsion equation, (2.?). 

All these equations assume small displacements and slopes, uniform 
density and distribution of stiffness, and the torsional equation is de- 
rived for a circular shaft. 


for Che „ plane bendln,: 

SA (21 

for the y~z plane bending: ■ - (EXL ^ VuJ 0 , 

w SA ^ (2 - 2 

where p is the density of the beam, A its cross sectional area, and 
(EDt x , (EI) y its (x-z) and (y-z) plane bending stiffnesses, respectively. 

Assuming separation of variables for u(z,t), one may write u(z,t) - 
r x( z >Px(t), and equation (2.1) can then be rewritten as: 


/A r x (2 ’ 3) 

This equation is valid if and only if both sides are equal to a constant: 


-0) 


Therefore P x + V 2 E^0(2.4), which integrates into 

“ cos C«rt + a), (2.5); where a is a phase angle. 

ih Yx = 0 r (4 LA* |J;0 

^ Wkre a*. M 

' x ~ (Eli 

this equation yields. 
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5? 


** = A '^h +& , C, fat fa + 4 Lntfa 
= 5 > UQ,t) = &>(aiiic(){A,fafe<8 l fa.} sCfafy^fa&jfa.n 

A si m ii ar reasoning enables us to find the solution of equation (2.2) in 
the following form: 


tyO= (2.8) 

Finally the z axis torsional bending is described by: 

(2 - 9 > % <; V&O 


f Zjz~ 

the modulus of rigidity of the beam. 

fuming ♦(«,« - 9(z) P 2 (t) znd substituting It into 
yields: 

P/ _ /^ ‘ 2ft) 


where G is 


equation (2.9) 


% - -f rw* 

5 J & 


( 2 . 10 ) 


awl 

* ?■ * * 


( 2 . 11 ) 


Therefore, 


- A ifa^z v 3j VZ. \l ) 


= enfant +£)f Ajfafo ] (i . I3) 
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B * Boundary Conditions (Z-X) and Natural Frequencies of Vibration 

The following relationships between shear, moment, and beam dis- 
placement are used in the boundary conditions 


_EI (3) 
L 3 u 


M — li v (2) 

x t2 


_ El 


Lr 

GI r 


M 


3 * 


3e 


(3) 


L 3 

- El u (2) 


and 


(2.17) 


Wbere^ V x * shear force in the x direction 
V y " " " y direction 

M x M y and M z the moment x,y, and z components, respectively. 
Ip is the beam polar moment of inertia. Let M,be the mass of the 
Shuttle while M r is the MS3 of the re fie c tor. The displacement in the 
direction of a point located at z * 0 is given by u(0,t)-Ay o $(0,t) 
and that in the y direction by v(0,t) + Ax o *(0,t) where Ax Q , Ay 0 are 
thecoordinates of the c.m. of the end body (Shuttle). 
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Now, an attempt will be made to cast the 10 equations describing 
the boundary conditions into the following matrix form: 

[ M ] (A} - 0 which has a non-trivial solution only when det [M] = 0. 

Since there is no offset at the Shuttle end, AXo » AYo = 0. Therefore 
B.C. (I) becomes 

rj’j = + Ms co* r,/ 

£ 3 U /fro 

Explicitly 

( £I /^)A i J-A, + C / j - ( v j 

B.C. (II) becomes 


— r A/ 

£3 J /£=ff 


At cO 1 r y/ 

A • / /£z O 



jj A J- A f + c, J 

=r- CC*A/rf_ 8/ 
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Equation IV* 

^ 4^ ¥ B^fafa + Cl “dfiz ¥ h fat fa / 


= ~ fan fa - ^ *2t^ - Ci^mtfa - ^ 

- A X L A s jtohfa - ^ £j 

Equation V' 



Equation VI' 

EZ/iff-B,*!),] = -a/]f + ZysfiM+Cjj 


Equation VII * 

§ 4 V- ^i/™A - % + Cl fat fa 4 bicfafaj - 

-<£*-{ Igrfr (Azcetft -fafa/4, tc t c**tfc 4 h ) 


— S 2 T &( Ayccofa C,Cto tfo + b/ fiutfajj 


Equation VITT 1 


'V”7 4, - S? cnA, + c, sJ,ty 3, y 4 tfo^ ; _ 

_ aytj ^ r fii (fy eafij 4 Ct cnJft 4 
- Zp 'Sfa(A,u>ft- %'>*& + C/ont ft y 
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Equation IX ' 


Ur /4 4 = - c# J }is 4 

Equation X ' 

ft* ( ^ ~ — - CD*£- fSjesyi, ) 

~i~ (- AXl £ A ' 2 +■ Bi Coife + Ci 1 •/■ 33 cnrf/%1 J 

•+ ^XM /&*/&/ -t B/ &F> fi, + C/ ft, + 3), c*rffo]J 
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3. GENERIC MODE EQUATIONS 


Consider an elemental mass, dm, of the body whose instantaneous position 
from the center of mass of the Shuttle is r. The equations of motion of 
dm can be written as 

a dm - L (q ) + fdm + edm ^ i) 

where a is the inertial acceleration of dm; f , the gravitational force per 
unit mass; e, the external force per unit mass; q,the elastic displacement 
of dm; and L, a linear operator which when applied to the small elastic 
displacement, q , yields the elastic forces acting on dm. 

The gravitational force per unit mass f, can be expressed as 
-+■ 

f * fo + Mr 

_ (3.2) 

where fo is the gravitational force per unit mass as the center of mass 
of the body considered and M q = matrix operator. 

In what follows, the generic mode equations will be derived based 
on a Newton-Euler formulation. The principal assumptions made in this 
development are: 1) within each component of the system, the mass 

and structural properties are uniformly distributed; 2) rhe material 
of each component is isotropic; 3 ) the system is deformed in such a 
manner that it experiences only small strains (within the linear range): 

4) elastic displacements are small as compared with the characteristic 
linear dimensions of the system; 5) the natural mode shapes of free 
vibrations of the system are known £ priori; 6) the system is nominally 
earth pointing; 7) the system is considered to be closed: no mass 
transfer across its boundaries. 


4<T 
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moving 


The vector equation (3.1) can be written in the frame 
with each body as: 

+ «»((&( r)]drn-L(^) , (f*e)d» (3.3) 

Note that r and r are the velocity and acceleration of dm as seen from 

the body fixed frame. The symbol m* refers to the inertial angular velocity 

of the body. The instantaneous position vector, r, of dm can be written 
as *— => ^-> 

(3.4) 

where r Q is the position vector of dm with respect to G, center of mass 

of the Shuttle, in the undeformed state; q* is the elastic displacement 
of dm. Hence 

^ •• .* 
f = o, and r* = a* 

1 1 (3.5) 

For small amplitude elastic displacements, one can write the elastic 

displacement, q, as a superposition of the various modal contributions 
according to 


q = Z 

where ~ (t) ( ^ + 6 *)'/<• 


(3.6) 

= modal amplitude 


and = r.t ^ 

V xk + ~ e* 

(3.7) 

The mode shape * (n) ( 2 ) is associated with the natural frequency, «o n , and 
satisfies the following conditions 
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(3.8) 


where M n is the generalized mass in the mode. 


L(f c ”) =. S coi q 1 ' 1 (3.,, 


( J w d m = o’ 

J M 

( f 0 * J dm - ~o 


(3.10) 


(3.11) 


This here assumes that the structural frequencies are much greater chan the 
1.745 hour/orbit * « 0 = 0.001 rad/s orbital angular velocity. This enables 
one to use, with a high degree of accuracy, the mode shape functions 
corresponding to non-rotating structures. The generic mode equation is 
obtained by taking the modal components of 'all internal, external and 
inertial forces acting on the system* i.e.* 

f r ^ ^ ^ 4 . -i 

l $ •[G<n» + r + <2 cox r t ODxr t co x (a3xr)J elm 

= l H $ C '[ L( fydm + T +? J dyn (3.i2, 

The various terms appearing in equation (3.12) can now be expanded as 
follows: 



(3.13) 


(3.14) 
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j M 1*3* Vdm r l £ | | ,( "!( ( ^q*Jcltr l 




( 3 . 15 ) 


16 ) 


$*! (<£*(&?)) dm = j ($'} S^Sxto)^ 


17 ) 


f ^ l LISJ dn\ - - cot 4* Mm. 

-'M dw 

fdnt = M? 0 d», 

4- f $ W * Mqdni < 


( 3 . 18 ) 


M 


<r I - 


€ dm - E^m 


( 3 . 20 ) 

where is the modal contribution of the external forces in the n ch mode. 
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/£dt(&L£ j A/\ 



Il r% = 


J// + In + 


In 


Iz-i tlrz + Minty 0 


hi 
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hi ■+ Irz 


“H = 


T/ 0 If 


0 
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0 
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I. MODELLING OF THE SCOLE CONFIGURATION 

• PARAMETRIC STUDY OF THE IN- PLANE SCOLE SYSTEM - 
FLOQUET STABILITY ANALYSIS 


THREE DIMENSIONAL FORMULATION OF THE SCOLE SYSTEM DYNAMICS 

• Rotational Equations of Motion 
Structural Analysis - Boundary Conditions 

* Generic Modal Equations 

n/ * WHAT WE CAN LEARN ABOUT THE OPEN LOOP SYSTEM? 


Consider SCOLE configuration without offset of the 
mast attachment to the reflector and without flexibility 

Consider SCOLE configuration without mast flexibility 
but with offset in-the direction of orbit (strawman) 

• Consider SCOLE configuration with offsets in two 
directions but neglecting mast flexibility 

Consider general SCOLE system dynamics 
♦ IMPLICATIONS FOR CONTROL STRATEGIES 


^ COL E f, /Vo A /g )' 

(£>[ 1 „ , *M% f Mri*+ 2 r]+ (pj /} . eQ, <jfi _ Tji t Ij} 

+1 T , - Irt + I r} J +c401^ 

- I 3 - Iu +Iri-Trt, t 3 (Ij~I z )-H>I^J '= O 

y Iil+ ( Pl T SS + lr *)+ & ( j / [lu + Isi-li l+ X r n./VI„Rl_MrRe 

+ Tri ~ Tr *] - $[Z-lu + hn- In - Mr£i] 

+ cDo (// (hz + 3 Zf ) — o 


9 9 

0/4 + In + Mri\ H.l%] + 3 a£ 9 ( I -I) _ 3ccfl< _ 0 

The. " 9 , pitch ’ Cpmti'cV 

Jkcoup&A Sum f&L dwo 

7 U 

,Oud /Zwt T,-I 3 <o and Zs+I**NrR ><? 
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F’jM.fkemont, Jkt Jaat ejaa^OT) /h£ ^oa: 
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fil 1 " <- M 4 * Tr,/ + (p fZ} y Afrxej- a 0[1, . Zt 

+ ?n +L, +I„ -ZtJ - a? 0(Zi* MrX£) _ oOty[T n 

+ Irj - T n - Kt t 3 (Zg _ fyJJ - 
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IMPLICATIONS FOR LINEAR CONTROL STRATEGIES 


After suppression of mast vibrations, linear system 
e q ns. have constant coefficients, control laws can be 
synthesized based on LQR techniaues. 

(A) For the special cases where the in-plane rotational dynamics 
separate from the out-of-plane dynamics, separate control 
laws can be generated for pitch and the roll -yaw systems. 

(B) When reflector offset results in coupling between the 
in-plane and out-of-plane systems, a bias momentum 
scheme could be considered so that the controllers serve 

to decouple the system via removal of the relevant coupling 
terms. Care should be taken so that saturation will not occur 

(C) Since the vibration frequencies of the mast are much greater 
than those of the gravity-gradient forced rigid rotational 
modes, actuators placed at strategic points on the mast 
could be used for quick removal of the vibrations without 
inducing substantial disturbances on the rigid modes. 

Once the mast deformations have been reduced to a specified 
level, the techniques described in (A) and/or (b) could 
than be utilized. 



II. CONTROL ISSUES: 

V 7 • CONTROL OF LARGE STRUCTURES WITH DELAYED INPUT IN 
THE CONTINUOUS TIME DOMAIN 

'Z ■ • CONTROL WITH DELAYED INPUT IN THE DISCRETE TIME DOMAIN 
^ ‘ C0NTR0L LAW DESIGN FOR SCOLE USING LQG/LTR TECHNIQUE 
• OPTIMAL TORQUE CONTROL FOR SCOLE SLEWING MANUEVERS 
Kinematical and Dynamical Eauations 
Optimal Control - Two Point Boundary Value Problem 
Estimation of Unknown Boundary Conditions 
* Numerical Results 

Discussion and Further Recommendations 
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I».B STABILITY ANALYSIS OF A SECOND OEDEE SYSTEM WITH DELAYED INPUT 

Th. vibration analysis of lar g . apac. structures la perform using 
« 0 NaL anal/aia ano sodal coordliates. transforning a coupled second ocder 

differential equations or partial differential equations Into n Oeooupleo 
second order differential equations of the fora 

*i + <*»l 2 Xi - fa 

... Cl) 

• •• 

where xj-i th modal coordinate 
^i™ 1 ^ natural frequency 

fi- influence of the actuators on the i th mode , and 
the control law of the form 

f i “ 2c i w i x i 

( 2 ) 

controls and stabilizes, the system (1) . The affect of delay In the 

control foro, was investigated with nuaerlcal emulation for the 

following numerical example. 1 

•• * 

x i * Sx^Ct-h) + 36x t -0 

TV (3) 

was Observed that for delay, h > 0.15, instability results,. 

The analytical verification of the above observation is obtained as 
follows^: 

The roots^of the characteristic equation 
G(s,h) - l p i(3) e -shi _ 0 

(4) 
i-0 

can be evaluated from the auxilary equation 
I p iCa)(l-Ts) 2 i (l«-Ts)2n-2i , 0 

(5) 
i-0 

where e'Jwh . fizW] 2 

l 1+ ^ T J (6) 
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( 7 ) 


( 8 ) 


Applying the above result to equation (3), the corresponding 
character 1 st inequation is given by: 

G(s,h) - l p ± ( 3 ) e"3hi 
i-0 

where P 0 (s) - s 2 + 36 
Pi ( 3 ) - 6S 
The auxilary equation is written as 

t2 S 4 + (2T ♦ 6T 2 ) S3 ♦ (1 ♦ 3 6t2-12T) S 2 
+ (72T * 6) S ♦ 36-0 

Using the Routh-Hurwitz criterion, Equation (9) has imaginary roots for 

T-0.0426 at <0-9.7. Using relation (6), h can be evaluated as: 
uih - ir/2 

( 10 ) 

or h s 0.16 

( 11 ) 

It 13 also brought to our .tt.nti.n3 that the above result oan be arrive. 

at without the approximation (6) for a aeoona or.er ayatem as follows, 

The characteristic aquation for ayatem (l) with the control law or 
the form 

f i - -2c i a> t x l (t-h) 
is written as 

s +. 2 _ g (13) 

To evaluate the minimum h for which equation (13) has unstable roots 
replace S by jo* as: 

-co 2 ♦ j2c i « 1 e"J“ h a) ♦ <u.2 . o 

1 (14) 

Using e"^* 1 - coswh-j sinuh, 

(15) 

Equation (14) can be written as: 


(xo 2 H- iCiCo^sincoh + cp ♦ J (ZCi^iCosaih) - 0 


(16) 
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Thus for equation (16) to be valid 
cosuih • o 

A ■ 

or uh -I (2P > i) 

P - 0,1,2,... 


( 17 ) 


and 

^“2C|W[(ii5iaj)Ji-oj2 m g 

the roots of Equation (18) are 
• • «i UiSinJh + /IT" Ci^} 

Takinff the Positive u and substituting into (17 ) 
. tt( 1+2P) 

Cl * 11 

2a) t { C^inaih+Zl+s^} 

Thus giving 

tynin " 0.1618 

for the numerical example (3). 


(18) 


(19) 


( 20 ) 

( 21 ) 


Thu, the example second order ayste. considered with the nature! 
period or oscillation or 1 second can not tolerate .ore than 0.16 seconds 
or delay without becoming unstable. Thus the general problem 0 r delay In 
control input must be carerully considered In the control system 
implementation of large space structures. 



the beginning. However, the delav in . 

delay in xnput in the discrete time domain 

can be relatively easily solved as shown below. 10 
The dynamic system described as: 


X(i + 1) * Z Aj X(i-J) ♦ I^BjU(i-j) 
J "° >1 


( 53 ) 


can be written as 

x(i+l) "1 

I 

x(i) 

x(i-m*l) 
u(i-l) - 
u( i ) 


u(i-2) 


. i 


u(i-i+i)| 


Z(i+1) 


A o A l‘**#A n B x B 2 ...b7 
1 0...0 0 0 0...0 
0...I 0 o 0...0 
0...0 0 0 0 . .00 
0...0 0 I o . .00 
0 o I. .00 

0 0 00 
L° °--*° 0 0 0..I0 


r 


0 

j 0 

; o 


X(i) 
x(i-l) 
x(i-m) | 
U(i-l) I 


U(i) 


i 


0(1 -2) , o 

! I 


a 

I 

L i«J 

Z(i) 


ORKifl'tAf, P 4 “ jo 

POOR QUALITY 


(54) 


(55) 


which can be written as: 

Z(i+i) . a Z(i) b U(i) 

Thus the augumented dynamic system (52) can be solved as a standard 

control problem. The only disadvantage is the increase in dimensionality 
o an already large dimensional problem. 
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a) Full State Feedback 



b) Observer Based Implementation 
- (sI-A ) _1 
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Important Properties of the Two Types of Implementations; 

1. The Closed Loop Transfer function matrices from command f to 

STATE X ARE IDENTICAL IN BOTH IMPLEMENTATION 

2. The Loop Transfer function matrices from Control Signal U' to 

Control Signed U (Loop Broken at xx) are identical in both 
Implementations 

3. The Loop Transfer function from Control Signal U" to Control 
U (Loops Broken at Point x) are generally Different. They 

ARE IDENTICAL IF THE OBSERVER DYNAMICS SATISFY: 

K II + C (SI-A]§ = BIC(SI-A) 1 b] for all S. 
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For Full State Feedback 

X = $ B U" 

For observer Based Implementation 
(I + $KC) X = $BU' -h^KC^BU" 

x = ( ♦ + KC) ( BU * + KC<t>BU") 

x = (I + <t»KC) (J> (BU' + KC<fr BU") 

= (I - ♦K ( I + C $K ) ^C) 4> (BU' + KC<|>BU" ) 
= 4>CB(C<t>B) -K ( I + C4>K) -1 ] C4>BU' 

+ ♦ [K-K (I + C4>K) -1 C<t>K] C0BU" 

= 4>[B(C«DB)~ 1 - K (I + C-HC) -1 ] C«j»BU' 

+ ♦K [I - (I + C«K)' 1 C'0K] C4»BU" 

= ♦[B(C0B)~ 1 - K (I + COK)” 1 ] C4>BU ' 

+ ♦ [K (I + C$tC) _1 ] C<1>BU" 



= [ I - A(I+BA)“ 1 B) 


An observer Adjustment 


Procedure : 


= 2 (q) C T R X 


AZ 


+ Z A + Q ( q ) _ 


T -1 

2 C l R l C 


Z = 0 


Q and R are treated as design Parameters 
fFor Kalman Filters 

, these are noise intensity 

matrices ] 


Q (?) = Qo + q 2 BVB T 
R = Rq 


For q=o 
For q -t- <D 


or 


K(q) is the nominal Kalman gain 

KRK T t 

— o— BVB 

q 

K ± ± , 

q - BV 2 (R 2 )~l 
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II. CONTROL ISSUES: 


• CONTROL OF LARGE STRUCTURES WITH DELAYED INPUT IN 

THE CONTINUOUS TIME DOMAIN 

• CONTROL WITH DELAYED INPUT IN THE DISCRETE TIME DOMAIN 

• CONTROL LAW DESIGN FOR SCOLE USING LQG/LTR TECHNIQUE 
v/ • OPTIMAL TORQUE CONTROL FOR SCOLE SLEWING MANUEVERS 

• Kinematical and Dynamical Equations 

• Optimal Control - Two Point Boundary Value Problem 

• Estimation of Unknown Boundary Conditions 

• Numerical Results 

• Discussion and Further Recommendations 
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OPTIMAL TORQUE CONTROL 
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Optical Torque Control 
for SCOLE Slewing Maneuvers 


PURPOSE : 

TO SLEW THE SCOLE FROM ONE ATTITUDE TO THE REQUIRED 
ATTITUDE, AND MINIMIZE AN INTEGRAL PERFORMANCE INDEX 
WHICH INVOLVES THE CONTROL TORQUES. 

CONTENTS : 

1. fCINEMATICAL AND DYNAMICAL EQUATIONS 

2. OPTIMAL CONTROL TWO-POINT BOUNDARY-VALUE PROBLEM 

(TP3VP) 

3. ESTIMATION OF UNKNOWN BOUNDARY CONDITIONS 

4. NUMERICAL RESULTS 

5. DISCUSSION AND FURTHER RECOMMENDATIONS 
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l . icic-U j.ii Dynoruc il Sguationa 

(.^ijiJ 3C0LC Configuration) 

1 = (1/2) T 1 

I V = - ,yj I v7 + \X 

v/:1 - ce 3 "’uier Parameter /actor : = fa . i , |T 

1 1 *0 - I *2 -^3 J 

w Angular VeLocity Vector ./ = (•,/. / f 


u Control .’orque 

3 -w, -Va, -w 3 

^ _ J i ’ Vs - v 2 

~' ; 3 J - 7 J 

, /z 7 


hi 

-I<z 

_I I3 | 

1 = -Ili 

taz 

■ Ii3 ! 

. - r l3 

~ l *3 

j 

Mere (Ref.l) 

11 = 1132533 

It 

v~* 

73 37 4 47 

1*. = -7555 , 

r 

ha - 

115232 

„ = 1535474 

' in* 

9533321 

iz = -13243 

' i <3 = 

155193 



1 53 - 71139 5 2 , 

7^3 = 52293 (Slug-ft*) 

I>3= 9545235 , 

I 2 3= 7 3933 (Kg-u i ) 



Transfer I to a iiajonal fora by an ortnoqonal natrix Z~*= 2 


3 . 93331 13 -3. JJ11151 3.3132393 
-J.3JS1684 3.9273353 3.3742533 
-'3 . 0132577 -3.3743342 3 . 9271252 


C T I C 


I 


3 A 


L >n 


wnera sub index , Tn , represents the principal axes system, 
I = 1133233 , I = 5935292 , I = 7137342 (Slug-ft" 1 ) 

Fro^ (2) , tine dynamical equation becomes 

C T I CC T w = -C T T CC T I CC T w * c T u 


or 


wner a 




"'m 


U = 2 u m , w = 

Similarily, we nave 

*n a ^/ 2 ) 

Eq.(3) can be written as 


(3) 


(4) 


(3) 


For simplicity, we drop subindex m in the following derivation. 
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/ytinai control 


Cost Function 


Two-?oint 3oundary-/ai ue Proulen (?P 3 \ 


J * n/uf* Si it = d/2) J ^ u T u Jt 

T ° -'to 


Tne Hamiltonian, H, for the system (4), (5) is 

H - (1/2) u T u + o T | + c t« 

&JSS.te,TJS* ta ‘' Principla ' e »« "<"««* conditions 
p = - O va i) === > p » (1/2) r? 


t = - (3 h/9 w} 


> r = U*/] r + ( 1 / 2 ) ( 3 ] p 


3 = PH/2>u} = = = > u = - f ' r 

plus (4) and (5), where o=[? D -> - f _ r t 

- l ^o - J , -° 2 P 3 1 » r = [r, r A r 3 ] T a 

tne costates corresponding to i and w, respectively 


[Jw] = 


J Z -'3 J 3 w 2 1 


T = ( I 3 “ 1 ^ ) / 


J, w 3 


fq] = 


j j ** ^ ■<, 


3 , -q, 


3 = ( it ~ I3 ) / 
J = (i 2 -i, )/ 


q 3 3 , -q 


After substitution of u from (3) into (5), we get 


— Jvviv — I r 


wnere 


Let 

2 *U, 

3 I 3 S 

//. W- 

1 4i 

V, 0 3 0 

3 ~o * 4 z 

- 3 3 r . 

c * r 3 f = (z 4 Zz) T 


V 


% 

1 # 2 a * 

fp. ?, 

P 3 P, r, Cj ] T 

Eqs. 

(4) ,(5) 

/ (7) ,(9) 

can be 

written as 





• 

2 3 

?(2) 



(1J) 


The boundary conditions 


2 i(t«) , z,(tj.) are known, 

2 a ( t e ) » 2 jt(tj.) ace unknown. 

Tnis is tne TP3VP. If we find tne unknown boundary values, 
z-(t,), then we can integrate (13) to get r, and fro.n (3) we 
obtain the control torque vector, u. 


7S 


- O - 
















3. -ouiTiation of Unknown 3ouniary Conditions 
3.1 special Casa of Slewing Motion 

Tha S "° t ' E cota,:as arbitrary axis Z fixed in both 

body axas 3 /steal and inertial space coordinate systen, i. a .. t ; le 

Sbler rotation. From the physical point of view, the rotation „ 
very ample, its rotation angle is small, and therefore nay 
.ones less energy (torgue). In view of oar C03t function, it 
is reasonable to think that the optimal slewing is near the 3 u- 
ler rotation. Considering the analytical solution about single 

2£i SCi?al aXi " -neneuver in Ref. 2, we define a rotation angle 
d(t) , aoout an arbitrary axis C, 


0( t) = 9 C + e 0 1 + ( 1/2 ) e 0 1 * + ( i/s ) o' t 3 


( 12 ) 


Foe the given boundary conditions 

e<J> **> ts,.g, (U) 

we have 

£= (5 ty/tJ)-<4 i/t f ) 

(?> -(12«f/tj) + (5 fc/tj) 

After substitution of 9 and S into <13., we can get z?b,.she 
initial guess of the costates at initial time t -t # . 


(14) 


3.2 Some Porperties of the Costates, 
Since q^q = i 


pT p » 3* = constant , but 3 2 * l 


we have 

3 is an unknown which is usually determined by iteration, thus 
[ Wj ] ===> 5 independent conditions 
t ?-i ^] T == = > 7 unknowns to be determined 
Fortunately, for the problem discussed in this paper, we can 
prove that 1 of the 4 unknowns o, can be arbitrarily selected. 
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Table I Slewing Data ani Boundary Values 



Fi&.l. UNE- 0F~ SI&HT EM04 


0 . 



Di3cu33ion and Partner Seco-xuiendations 


ConsiJar the Distribution of u on the Snuttle and 
Ref lec tor . 

riTie-Optinal Slewing, (Rigid Configuration), 

Cost Function 



Solve the TP3VP . by Snooting Methods 
Include the Flexibility in the Problens. 


tn 


^ f 'Jd I2 ^ 1 ^3 W3 Aj Ajj 

( lx 14 + 2n ] 
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O 


a 

w i 
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“2 


3 :Jo * °f flexible appendage nodes included 
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c= 


COS9 2 COS0 3 


-cose 2 sin 03 


sin© 


(sin 9l sine2Cose3 (-sin 9 jsin 92 sin 03 -sin 9j cos 9 2 
+sin 9 cose ) +cos9 COS0 } 

(-0036:3109200303 (0030:3100251003 CO30,COS 92 

- + sine 3Sinei) -003935100:) 


a: 


Srvlrto^’flSdTn^J :h h % deX , tral set of ^hogooa; 
the rotation Jf ^ eV t ^S#- Xed I^ me ' the " l 1 
rotation of k. 2 the rotatlon of j and e 3 is^he 

transformed from 1 the 7ne°rtfal ?» the f orbit er can be 
frame for the body-three angles as the body - fixed 

U) = M^g 

~ ( 2 ) 

can be h gi t V e t n al as k [ 4 ] t1C expression of the system 


T = T o + T, - T 2 


(3] 


where T 0 i s the kinetic enerav of f-ho 4 . 4 .*, 
given as e or the shuttle and i 


T 0 = 1/2 m: V T V + 1/2 iu T i x u 


(4] 


The kinetic energy of the flexible beam is T , and t 


95 - 


V# 



l/2mV^V + 1/2 aJj (j- mV J Cw+ 1/2 d T d dm 
+ V o T/d dm + wj/^ad dm + 1/2 [u x 'u Ju ^ ]dl 


u ' 

y 

u/ 



f 


T, = l/ZmVjo 1 V_+ l/2u T Ju- raVjj T cu+ q. 2 +Vg T £ 

t* „ n # o n # o 

+ “. S.+ l/4p C E p q f + ] 

i=l 5i l i=i bl 1 

where 


u x- 1 £ 1 *xl( s )‘«l( t ) 


n 


u y = 1 ^ 1 +yi (s)q i (t) 


n 


and 


u -i» ■ 1 £ 1 **i (s)q i (t) 


P = /* xi (s) ds 
li o 

L 

P 2i = / <*> yi ( s ) ^ 


0 


P 3i = / S ^xi (s) dS 
X 

p = J S(j> . (s) ds 

4i 0 yi 

P„ = / ( S *i ) 2dS 

5i o X1 

L 

p 6i “ / ) 2ds 


( 6 ) 


(7) 


( 8 ) 
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and 



n . n 


i=l P li q i 

« (t) = 

i-1 P 2i q i 

■ 

0 

L (t) = 

i-1 P 4i q i 


n 

p 

ill P 3i q i 


0 


( 9 ) 


(10) 


The kinetic energy T 2 , of the tip mass (the reflector) 

T 2 = 1/2 m 2^o ^0 ' m 2 ¥o!( L )“ + m 2 ^ ^ L ) 

- 1/2 (L)a(L)u + m_ «Ta(L)d_(L) 

+ 1/2 m 2 d_ T (L) d_(L) + 1/2 n T I 2 S. 


( 11 ) 


where 


ft = 0) + 


u x '(«■) 

“y 'W 

L% W 


(12) 
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T = 1/2 ra 2 v/Vj, - m 2 vj afLJu + m 2 vj d_(|_) 


- l/2nfew T a(L)aJL)c i + m 2ii T a(L)d^L) 


n n 


+ 1/2 ^2 C^i j!fxi (L) *xj (L) ^.^j + i I 1 ^ ( i L) *yj (L )^i'’^ 
+ 1/2_P_ T I 2 _P_+ 1/2o) T I 2 « J ’ (13) 


where 
n 


n i 


n 


(14) 


-“^-l x 1 (L,q i(t) i (t) 1 (t) ^ 

-instituting Tq» T 1 and \ from the foregoing equations 

^o equation (3), the total kinetic energy expression 
can be written as 


T = 1/2 m 0 V_ T V_ + uJ HV + I/2J I 0 v_ T A j g_ 


+ 21 A & + 1/2 Q A j q 


(15) 


where 


?? 




(26) 


%V_~ Hi ♦ A lf = + 

where the nonlinear terra Niis given as 

U 1 = C ( m o_y ■ H + A 1 q ) ( 27) 

* (in n V - H u+ A x _q ) 

Similarly, using equation (2) and the chain rule in the 
Lagrange's equations, the rotational equations are 
obtained as 

HV + I„ u+ A 2 q = G(t) + N 2 (28) 

* 

where jGft) is the net moment about the mass center of 
the orbiter and is given as 


G = G +(r + a)xF 
— - o — — —2 


(29) 

and the nonlinear term is given in terms of 

t ransf ormat i ons M and C, and ^ and .The vibration 
equations of the beam can be obtained by again using 
Lagrange s equations and the potential energy function 


U = 1/2 q '« q 

'■'Here the stiffness matrix K is given as 


(30) 


K = 


3r, 4 


IL L J j 

L 3 


(31) 


he vibration equations are 
. T* T- 

A i_V + Aj u + A 3 q = -Kq 


(32) 
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Io 0L + A 2 q = G(t ) + H 2 ( “ ) 

T. ~ ” 

A 2 w + A 3 q = -Kq 


Equation (33) can be rewritten as 

i = i 0 ' 1 [ g + 1 2 ( “) - A 2 q ] 

The first three Euler parameters are defined as 


£2 = V 4>/2 


= cos ^ /2 

= 1/2 ( £4 £+ £ x £ ) 


de_4 

dt 


= - 1/2 £•£ 

= 2 ( £4^ - £4£ 

d£ , 

= hT- » h. ( £, to) 


de 

e x ~ ) 
dt 
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VECTOR DIFFERENTIAL EQUATION FOR DISTRIBUTED STRUCTURES 
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Assume that the shuttle is fixed and identify the mast parameter. Use 
perturbation technique in conjunction with frequency response. 
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PARAMETER IDENTIFICATION (CONT'D) 
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PARAMETER IDENTIFICATION (CONT' D) 

If the measured output is not the whole state, use Kalman filter to 
estimate the state. First estimation is based on the postulated model 
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On Incorporating Damping and Gravity Effects in Models of 
Structural Dynamics of the SCOLE Configurat ion 


ABSTRACT 


The damping for structural dynamics models of 
flexible spacecraft is usually ignored and then 
added after modal frequencies and mode shapes 
are calculated It is common practice to assume 
the same damping ratio for all modes, although it 
is known that damping due to bending and that 
due to torsion arc different Mass effects on 
damping are sometimes ignored 

It is i-he purpose of this paper to examine two 
ways of including damping m the modeling process 
from its onset First, the partial derivative equations 
of motion are analyzed for a pinned -pinned beam 
with damping. The end conditions are altered to 
handle bodies with mass and inertia for the SCOLE 
configuration Second, a massless beam approxima 
tion is used for the modes with low frequencies and 
a clamped -clamped system is used to approximate 
the modes for arbitrarily high frequency The 
model is then modified to include gravity effects 
and is compared with experimental results 


0 


• Introductory Remarks 

• SCOLE Configuration 

• Partial Differential Equations 

• Pinned-Pinned System with Damping 

• Free-Free System with End Bodies & Damping 

• Massless Beam Approximation 

• Gravity Effects 

• Comparison of Model Frequencies 

• Concluding Remarks 
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Shuttle (and Reflector) Body 

*1 = (F 1 * f'l.Beam)/™, 

t, t = -ffij, 7 

Roll (and Fitch) Beam Bending 

PA 0^2°~ CI 0^T° + EI d 

dt ds^dt ^ ds^ 


Yaw Beam Torsion 


D f 0.n 6(s ' s n> 
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+ 8 ° n 5s 
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Beam Elongation 
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Damping Considerations 


• The Classical Damping, ^ — Yields Excessive 

^ _I 4 


ds dt 

Excessive Damping at Higher Mode Numbers 

• The Term, 2 is Consistent with experimental 

ds dt 

Data. 

• The Practice of Post- Analysis Addition of 

Damping Ignores Effects of Mass. Stress Type 

• Damping Must be Included from the Start. 
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Distributed Parameter Model of 
SCOLE with "Proportional Damping' 


# Start with Pinned-Pinned Beam with Damping 

• Add Bodies with Inertia at Ends 


• Model Acceleration of Frame as Inertial Loading 

• Extend in Three Dimensions to 

SCOLE Conf iguration. 

• Yields Infinite-Order, Modal, State Equations. 
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Distributed Parameter System 
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|^tassles^Be^imModeI| 

• Exact Static Deflection 

• Approximates Low-Frequency Modes 

• Nonlinear Kinematics 

• Linearized State Space, Modal Model 

• Classical DampinglWorking Proportional) 

• Extended to n-Body Network 


131 




^Jravit^^ffectsl 

• Assume Cubic Deflection of Beam 

• Express Potential Energy due to the 

Raising of End Body 

• Relate to Stiffness Matrices of the 

Massless Beam Model 

• Incorporate Gravity Effects in the 

Stiffness Matrices 

• Gravity Effects Larger than Structural 

Stiffness 
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Equations of Motion 
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a32 = 5 1 I_F U i l +I kl = - a 36 

a 34 = mjt F u ] = a 38 

a 52 - r,'l-M u r 4 . r k r; l- u r ]• - a 56 

a 54 " 4 l M u'5 F ul" ~ a 58 



I - Moment of Inertia 
m- Mass 

r - Coordinates of attach point 
r - Cross product operator, rx 

Mu.M^,F u> 5: - Stiffness Matrices 
1 - Denotes the Shuttle body 
4 - Denotes the reflector body 
u,^ - Beam deflection and slope 
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c 


c 


c 


c 


program Tuoeoo 


ORIGINAL PAGE IS 
Of- POOR QUALITY 


REAL I 1, 14, HIM, KIN,m,n4,nU,MRNG,L,nHSS1,MRSS4 
DIMENSION IK13), 1 HN<13), I4<13), I4IN< 13>,A(580>,RflTl<?> 
♦ROT 1T< 13 ), RRT4<? >, RRT4T< 13 >, FU< 13 >, FAN0< 13 >, MU< 13 >, HANG< 13 ) 

*OUM(580),OIJ1<588),EREflL<38)*EinAG<38>,EUEC(580>,DUCK580) ^ 

- . ... DEFINE INERTIA HRTfilCES. .... 

CflLL SET(A,24,24> 

CALL SET<I 1,3,3) 

• 1 <5 >*905443. 


I 1<7)=- 145393. 

11 <9 >=6789 100. 

1 1< 11 >—145393. 

I K 13 >=7886601. 

CflLL SPITO 1,3H ID 
CflLL MflKE(DUH, M) 
CflLL INUR< 1 1, 1 11N> 
CflLL MAKEd 1,0UN> 
CflLL SET< 14,3,3) 

I 4(5 >=4969 
I 4(9 >*4969. 

I4( 13)=9938. 

CflLL SPIT< I4.3H I4> 
CflLL HflKE(0Un, I4> 
CflLL INUR( 14, « 4 1 N > 
CflLL MflKE(l4,DUM> 


CflLL SET(RflT1,3, 1> 

CflLL TILDfl(RAT1,RflT1T> 
CflLL SET(RflT4,3, 1) 
RAT4<5 >=- 18 75 
RflT4<6>=32.5 
CflLL TILDfl(RflT4,RAT4T ) 
MflSS4= 12 . 42 


DEFINE ATTACH POINT UECTOR, MATRIX 


tW ^ T ° REFLECT ° R BOOV 

CALL AOO(AD,RftT4T,-1 . , RflT4T 0UM> 

CflLL SPIT (DUN, 4H DUN) 

CflLL NULT(DUN,DUH,OUN) 


C°LL SP!T(QUM.4H DUN) 

CflLL flDOd., 14,-12. 42,DUN,I4> 
flO= 5*12. 42/( 12 42+ 5*12 42 > 
CflLL RDD<RD,RRT4T,-1.,RRT4T,DUH) 
CflLL f1ULT(DUH,0UN,DUN) 

CALL flOOd , 1476 21, DUN, 14) 

CflLL SP.'T(I4,5H i4NU> 

CflLL INUR< 14, KIN) 

NflSS 1=6366 . 46+ 09556* 130. ( 2 . 
NflSS4=HflSS4+ 09556* 130/2 . 


El =40000000. 
GJ=40000000. 
Efl= 100000000. 
L=130 


BERN SECTION CHARACTERISTICS 


CflLL SET(FU,3,3) 

FU<5 >=- 12 . *£ I /(L*L*L ) 
FU(9 >— 12 . *E I /<L*L*L > 
FU(13)=-Efl/L 


CflLL SET(FflNG,3,3> 
FANG ( 6 )=6 *E I /< L*L ) 
FflHG(8)=FflNG<6) 


SET UP FORCE /OEFLECT I ON MATRIX 


SET UP FORCE/SLOPE ANGLE MATRIX 


C 


CflLL MAKE<NU,FAf1G) 


SET UP MOMENT /DEFLECT I ON MATRIX 


c 


SET UP nOnENT /SLOPE ANGLE I1ATRIX 


♦ 


C 


c 


CflLL SET<t1flNG,3,3) 

HANG<5)— 4.*EI/L 
nHNG(9)=-4 . *E I /L 
HflNGC 13>=-GJ/L 
CflLL SPIT<FU,3H FU) 

CflLL SPIT<FANG,5H FANG) 

CflLL SPIT<HU,3H HU) 

CflLL SPIT<nflMG,3H flflNO) 

CflLCULflTE ELETEMTS 

CflLL HULKRflT 1T,FU,OUH) 

CflLL MULT<DUH,RRT1T,OUN) 

CflLL HULKRflT IT, FANG, OUH) 

CflLL flOOd., Dlt1,-1., OUH, OUH) 

Cfla flOOd., NANG, 1.,0UN,0UN> 

CflLL HULKHU,RRT1T,OUH) 

CflLL AOOd1.,OUH, 1.,OUH,OUH> 

CflLL HULK1 1 1H,OUN,OUM) 

CflLL INSERT( t,4,DUH,fl) 

CflLL flOOd 1.,DUH,0., OUH, OU1) 

CflLL INSERTS 1, 16,0UH,A) 

CflLL HULKRflT IT, FU,DUH> 

CflLL flOOd, MU, 1, OUH, OUN) 

CflLL HULK 111 N,OUH,OUH> 

CflLL IMSERTd, 10,OUH,fl) 

CflLL flOO<- 1 ,DUH,0. ,OUH,OUH) 

CflLL 1 NSERT < 1 , 22, OUH, fl ) 

CflLL I0ENT(DUH,3) 

CflLL INSERTC4, 1,OUH,A> 

CflLL I HSERT < 10, ?, DUH, fl ) 

CflLL INSERK 16, 13,DUH,fl) 

CflLL INSERK22, 19,0UH,fl) 

CflLL MULT <FU, RflT IT, DUH ) 

CflLL flM3d1 . ,OUH, 1 . , FANG, DUN) 
flO-l./HflSSI 

CflLL flDD<flO,0UM,0. ,DUM,OUH) 

CflLL INSERK7,4,DUH,fl) 

CflLL flOOd 1, OUH, 0., OUH, OUH) 

CflLL INSERK?, 16,DUH,fl) 

CflLL AOO<flO,FU,0.,FU,OUH) 

CflLL INSERK7, 10,DUH,fl) 

CflLL ftOOdl . ,OUH,0. , DUH, DUH) 

CflLL INSERT<7,22,0UM,fl) 


CflLL MULKRfiT4T,FU,DUH) 

CflLL MULT<DUM,RftT4T,OUH) 

CflLL MULKRAT4T, FANG, OUH) 
CflLL flOOd., OUH, -1., DUH, OUN) 
CflLL flOO< 1 . ,DUH, 1 . ,MflHG,DUN> 
CflLL HULKMU,RflT4T,0UH) 

CflLL flOOd 1 . , DUH, 1 . , OUN, DUN ) 
CRLL HULK 1 4 IN, DUN, DUM) 

CflLL INSERK 13, 16, OUH, fl) 

CflLL flOOd 1 , DUH, 0. , OUH, DUH) 
CflLL INSERK 13, 4, DUH, fl) 

CflLL MULKRflT4T,FU,0UN> 

CALL flOOd , DUN, )., HU, DUN) 
CflLL HULK I4IN,DUN,DUH) 

CflLL INSERK 13, 22, OUH, fl) 

CflLL flOOd 1., DUH, 0., OUH, DUH) 
CflLL INSERK 13, 10, DUH, fl) 

CflLL HULT <FU, RAT4T, DUH > 

CflLL flDO<-1 . ,DUH, 1. , FANG, DUN) 
AO= 1 . /HASS4 

CflLL flOO<AO,OUN,0.,DUH,OUM> 
CflU INSERT< 19, 16,OUH,fl) 


IN "fl" MATRIX 


ORIGINAL PAGE fS 
0F POOR QUALITY 
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£*-l Aoo<-i.,oun,0.,oun,oun> 

CflLL IMSERT< 19,4,0UM,A) 

CALL AOO<AO,Rj,0.,nj,OUH) 

CALL IMSERTC 19,22,DUTI,A) 

CALL AOO<-1.,OUH,0.,OUH,OUH> 

CALL IMSERT<19, 10,OUM,A> 

OLL E i GOKA ‘ EREAL ’ E I MAG,' EUKJKA^ SHflPES 

CALL SPIT<EREAL,3H REAL) 

CALL SPIT<EIMAG,5H I HAG) 

123 F0ANAT< 1 10,E15.6) 

C 66’ 10 1*4^580 PR,NT NCW ~ ZER0 EL£nENTS ^ "B" MATRIX . 

IF<A<l)*^-.0000000000 1)11 ,, 12 
12 PRIMT 123,I,A<|> ' ' 

11 CONTINUE 
10 CONTINUE 
END 

— EO I /TOP— 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


M 1 2 3 

1 9054E+06 00006+00-. 14546+06 

2 . 00006+00 67896+07 .00006+00 

3 1454E+96 oooe o ee 7087E+07 

•43418836+28 

>4 1 2 3 

1 .49696+04 00006+00 .00006+00 

2 . 00006+00 49 ME +04 . 00 0 06*90 

3 . 00806+00 . 00006*00 .90386*04 
. 2453788E+ 12 

0U1 1 2 3 

1 .00006*80 - 0OO OOO O .10836*02 

2 .0000E+00 . O 000 E t 06- .02306*01 

3 .1083E+02 .62506*01 .00086*80 

OUH 1 2 3 

1 - . 1 1746+03- . 677 1E+02 . 00006+00 

2 - 677 1E+02-. 39066*02 .00ml£+00 

3 0000E+00 . 0000E+00- . 15646+03 

I4NU 1 2 3 

1 . 9342E+04 .25236+04 .00006*00 

2 .25236+04 6424E+04 .00006*00 

3 .00006+00 .00006+00 . 1577E+05 
84589596+12 

FU 1 2 3 

1 -.21856+03 .00006+00 .00006+00 

2 00006+00- . 2 185E+03 .00806+00 

3 .00006+00 .00006+00-. 76926+06 

FANG 1 2 3 

1 0000E+00 .14206+05 .00006+00 

2 .14206+05 .00006+08 00006+08 

3 0O00E+80 00006+00 0000E+00 

HU 1 2 3 

1 00006+00 .14206+05 .00006*0 0 

2 1420E+05 . Q068 E +08 6868 £ *00 

3 0000E+00 . 0000E+00 . 3000E+00 

MflNG 1 2 3 

1 -. 1231E+07 O000E+00 .00006+00 

2 . :o30E+00-. 1231E+07 0000E+00 

3 0000E+00 00006+00- 30776+06 

°ERL 1 

1 -. 1344E-04 

2 - . 9690E-06 

3 -.236 IE-08 

4 - 621 IE- 14 

3 S211E-14 

6 28616-08 

7 .96906-06 

8 . 1344E-04 

9 .0000E+00 

10 0000E+00 

11 .00006+80 
12 .00006+00 

13 .00006+00 

14 .00006+80 

15 . 00006+00 

16 .00006+00 

17 .08006+00 

18 .00006+00 

19 .0000E+00 

20 .00006+00 
21 00006+00 
22 00006+00 
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ORIGINAL PAGE IS 
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Static 


Number of Modes 



ctioti 


Error in Deflection 

39 % " 

24 

17 

13 
1 1 
9 
8 
7 
6 


6 
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Comparison of Modal Frequencies 


Mode 

P.D.E. 

Finite El. 

Lumped Mass 
*Clampcd 

1 

.278 

.276 

.258 

2 

.314 

.301 

.370 

3 

.812 

.810 

.926 

4 

1.18 

1.18 

1 79 

3 

2.05 

2.03 

2.57 

6 

4.76 

4.77 

4.28* 

7 

5.51 

5.52 

4.28* 

8 

12.3 

12.4 

1 1.89* 
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SCOLE Experiment 

Modal Characteristics 


Small 

Amplitude Amplitude 



O O Q 

5th Mode 


□ □ □ 


4th Mode 


Proportional Damping 
(Constant Damping Ratio) 


3rd Mode 
<£ 0 ~ 


-.02 


Real Part, tr 


2nd Mode 


list Mode 



[^oncludin^^en^rl^^J 

• An Infinite-Order State Space Model 

was Developed which Incorporates 
' Proportional" Damping. 

• A Lumped Mass Model of SCOLE was 

Developed which Includes Gravity 
Effects and Classical Damping. 
Extended to n-Body Modeling. 

• Modal Frequencies are Compared r«r 
the SCOLE using Different Methods. 

• Items Remain to be Addressed Before 
SCOLE Modeling is Complete. 
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MODEL DESCRIPTION 


DEFINE ELEMENT PROPERTIES : 

- MATERIAL CONSTANTS 
MODULUS OF ELASTICITY, 

POISSON’S RATIO, 

DENSITY 

- MAST, REFLECTOR, RIGID LINKS AS BEAM 

ELEMENTS 

- CABLE AS BAR ELEMENT (AXIAL STIFFNESS 

ONLY) 

- SHUTTLE AS VERY STIFF BEAM (ASSUME 

RIGID) 

JOINT LOCATIONS AND CONNECTIONS : 

- 44 JOINTS TOTAL, 7 FOR CABLE, 12 FOR 
MAST, 6 FOR REFLECTOR AND REST FOR 
RIGID MASSES 


IS) 




MODEL DESCRIPTION (Continued) 


TWO BOUNDARY CONDITIONS MODELLED: 

CASE 1— SUSPENDED (6 DOF FOR ALL JOINTS 
EXCEPT TOP OF CABLE) 

CASE 2- CANTILEVERED CABLE, SHUTTLE 
PLATFORM FIXED IN ALL DOF 


INCLUDE RIGID MASSES AND CONNECTIONS 

- ACTUATORS 

- SENSORS 

- SHUTTLE PLATFORM AND COMPONENTS 


CALCULATIONS : 

- STIFFNESS AND MASS MATRICES 

- INITIAL STRESSES (DUE TO GRAVITY 

LOADING) 

- STATIC DISPLACEMENTS AND REACTIONS 

- EIGENSOLUTIONS - FREQUENCIES AND 

MODE SHAPES 


FREQUENCY DATA FOR CANTILEVERED CASE (FIG 1,2) 


MODE 

FREQ (HZ) 
EAL LAB 

DELTAS ‘ 

EAL/LAB RATIO 

i 

0.443 

0.44 

0.7 

1.01 

2 

0.447 

0.44 

1.6 

1.02 

3 

1.504 

1.54 

2.3 

0.98 

4 

2.913 

3.00 

3.0 

0.97 

5 

4.345 

4.36 

0.3 

0.99 

6 

6.821 

3.08 

121.5 

2.21 



FREQUENCY 

DATA FOR 

SUSPENDED CASE 

(FIG 3,4) 


FREQ (HZ) 



MODE 

EAL 

LAB 

delta-:: 

EAL/LAB RATIO 

6 

0.566 

0.55 

2.9 

1.03 

7 

0.638 

0.65 

1.8 

0.98 

8 

1.514 

1.62 

6.5 

0.93 

9 

2.940 

3.10 

5.0 

0.95 


RATIO OF SUSPENDED TO CANTILEVERED FREQUENCIES (FIG 5,6) 


MODE* EAL LAB 

1 1.28 1.25 

2 1.43 1.48 

3 1.01 1.05 

4 1.01 1.03 


* NOTE: SUSPENDED MODES 6-9 CORRESPOND TO CANTILEVERED 
MODES 1-4 















COMMENTS 


- EAL, LAB DATA IN GOOD AGREEMENT 

- HIGHER MODES TEND TO HAVE SLIGHTLY 
LARGER DIFFERENCES BETWEEN EAL ic 
LAB RESULTS 

- FOR HIGHER MODES, FREQUENCIES OF 
THE SUSPENDED AND CANTILEVERED 
CASES ARE SIMILIAR; THE MODE 
SHAPES ARE ALSO CLOSE 


IGO 


CONCLUSIONS 

- EAL, LAB FREQUENCY DATA MATCH WELL 

- NEED TO GET MORE ACCURATE MEASURE- 
MENTS FROM LAB, AND WITH MORE 
MODES FOR BETTER COMPARISON 
COMPUTER MODEL Si LAB 

- FOR HIGHER MODES, THE CANTILEVERED 
CONDITION MAY BE SUBSTITUTED FOR 
THE SUSPENDED, THUS REDUCING THE 
NUMBER OF NODES AND DOF’S IN 
COMPUTATION 
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MODEL REFERENCE CONTROL 
OF DISTRIBUTED PARAMETER SY: 
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MODEL REFERENCE CONTROL OF LUMPED LINEAR SYSTEMS 



INTRODUCTION 
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LUMPED MODEL 
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DEFINE IDEAL STATE AND CONTROL 



ASSUME 
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APPLY 
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SPECIAL CASE (PMF) 
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SCOLE APPLICATION OF LUMPED MODEL FOLLOWING 
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SPECIAL CASES: 8 outputs 
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COMMAND 


(cont. 
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CA SE II (a ): Same objective as in CASE I (a) 



179 


In both cases: 4 states 
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MODEL (cont.) 
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SYSTEM I .C. ' s : RESULTS 
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- 0.10 - 0.02 0.06 0.14 0.22 0.30 


PL.BaJT STATE : C&5S- 



\% 


MODEL REFERENCE CONTROL OF DISTRIBUTED SYSTEMS 
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CONTROL PROBLEM FORMULATION 
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SOLUTION TO DPS MRC PROBLEM 


A 






}99 



U) q(t) 



200 



2 01 



FOR ILLUSTRATIVE PURPOSES WE WILL CONTROL 
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REFERENCE MODEL 
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CASE 2: Simple-free (shuttle mass infinity) 



CASE 1 Free-Free 



2 06 



SOLVING THE EQUATIONS OF MOTION 
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Purpose: To damp out the structural vibrations within ten 

SECONDS WITHOUT VIOLATING THE CONTROL MAGNITUDE 
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MODE » n k W n (rad/sec) j f (HZ) 
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A MODEL FOLLOWING PROCEDURE USING CGT THEORY WAS 
DEVELOPED FOR APPLICATION TO DPS SYSTEMS. 
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PLANNED ACTIVITIES 
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STATEMENT OF THE REGULATOR PROBLEM 
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OBJECTIVES OF THIS ANALYSIS 
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.02 degrees LOS pointing 



PROCEDURES OF THIS ANALYSIS 
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nruster moments on shuttle 
Thruster forces on reflector 
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DisCDS Sifu-ATiow: Rigid Bodies Connected by Hinges 
ftoss aid Mow of Inertia Provided by User for each Body 
Location of Hinses and Sensors 
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COMPONENTS OF THE ANALYSIS 
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control algebraic Ricatti equation (CARE) 

No noise or time delay in sensors or actuators 


LQ REGULATOR FOR FLEXIBLE BEAM 
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LQ REGULATOR (CONTINUED) 

♦ Objectives in minimizing cost functional: 

1 . Maximize reaulator performance 
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u(t) is recalculated each time step in DISCOS 
from current value of x(t) and offline, 
predetermined values of R, B, and P. 


TABLE I: MODES OF THE SYSTEM AS COMPUTED BY NASTRAN 


MODE NUMBER MODE TYPE 


ANGULAR FREQ 


FR eq. 


IN HZ 
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11 

12 
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PITCH 

ROLL 

YAW 

ROLL 

PITCH 

ROLL 

PITCH 

ROLL 

PITCH 

COMPRESSION 

ROLL 

PITCH 


1. 746 
1 . 969 
5.105 
7.410 
12.848 
29.459 
34.263 
74.670 
78.883 
106.281 
142.467 
145.618 


0.278 
0.313 
0 . 182 
1.179 
2 . 045 
4.689 
5.453 
11.884 
12.555 
16.915 
22.674 
23.176 
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INPUT MATRIX B 
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ratio of actuator influence to thruster influence 


TABLE 2: ACTUATOR LOCATIONS FOP. MAXIMUM CONTROLLABILITY 



Controllability surface, 
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excited bv minimum-time RAPID slewing 
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INTRODUCTION 

HIGHLIGHTS OF NUMERICAL RESULTS 

MINIMUM-TIME RAPID LOS POINTING SLEW FOR SCOLE 

ADAPTATION OF LOS ERROR EXPRESSION 

CONCEPT OF "MODAL DASHPOTS" 

MODAL-DASHPOT VIBRATION CONTROLLERS 
DESIGN AND SIMULATION RESULTS 

CONCEPT OF "MODAL SPRINGS" 

MODAL-SPRING U I BRAT I ON CONTROLLERS 
DESIGN AND SIMULATION RESULTS 

COMBINED USE OF MODAL DASHPIOTS AND SPRINGS 
MORE DESIGN AND SIMULATION RESULTS 

CONCLUSIONS 
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kL.N E- OF- SIGHT ERROR - - GENERAL VECTOR 


EXPRESSION 


0 FR = ray of emission = R R -R f 

FIR = REFLECTED ray = LOS VECTOR = R L0S 

° BT = FR + 2 RR* SINCE 

FR +RF“ = FT ; FR H-RR* = £R» ; FF* - 2FR* 

0 RR' = CRF . R a )R A = - C FR . Ra>Ra 


IN UN-NORMALIZED FORM: 


K LOS RF J - Rr Rp ~2£(R R -R f ) . R A J R A 
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TRANSFORMING TO INERTIAL FRAME 
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<T lRLOS>x] + [CTtRLosJyj' 




MORE ON LOS ERROR EXPRESSION 
INCLUSION O F MOST BENDING ftND TORSION 
Rr -Rp = R T - Tj T 4 R b -R p 
WHERE 


R*p 


BENDx 

BENDy 

“ \/l30* 1 BENDx 2 -BENDy 2 


Rb 


18 . 75 
-32 . 5 

0 


Rp 


BENDx = u (4) -u ( 1 ) 

^ X 


3 . 75 
0 
0 

BENDy =Uy(4) -u y (l) 


LOS x = - C T ^ R, \ - _ t' 

x i "Los J y - T 


lry R LOS 


2T 


4yz T 4xz’ I+ 2 T 4 yz' ~ 2T 


T 

4yz 4zz 


+ T 4r*y «B 


:J T 1 |[ R T R fJ 


LOS y = CT t R los ) x = T 


R 


lrx LOS 


I - 1 2T 4 xz' 2 T 4xz T 4yz' 2 T 4xz T 4zz] T 1 [ R ' 


R F ] 


+ ^4px T b 
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DYNAMICS: 


force 


fl 


d 2 


dx 


at 2 - - at 

(TORQUE) ACTUATORS 
f = Bp u 


+ K x = f 

AND VELOCITY SENSORS: 


MOW1AL modal representation x = 5 „ . 

— 2 *1 dR 

dtz A dt + «* >1 = §T Bp U 

WHERE 

ft 2 = DI AGJjwj 2 J = 5T k g 

A = d 5 


CONTROL LOW FOR CONSTANT- GO IN 

VELOCITY- OUTPUT FEEDBACK: 


FULL-ORDER CLOSED-LOOP SYSTEM EQUATION: 

a 2 1~| 

+ CA + « T Bf G Cy 5 ) + ft* 

V .. dt “ M 

A* 


0 
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MODAL -DASHPOT APPROACH 
DESIGN TO ACHIEUE I HD K PENDENT 

DAMPING augmentation for each mode in a 
reduced-order model 

LCT «i BE damping ratio desired of modeled MODE i 

SET 5 M T Bp G Cy = DIAGf~2CiWi"| 

v — nf 

THEN SOLUE FOR FEEDBACK GAIN MATRIX G, 

G = <5 M T B F ) f DIAGpCj^i] (Cy S|i>* 

USING THE PSEUDO- INVERSES ( )t DEFINED AS FOLLOWS 

<5, 1 T B F )f = <5 n T Bp)T [(5 M t Bp) (5 M t b F )T]-» 

* 

(c v S M>* = SS M )T (Cv Sr1 )J-« (Cv afl)T 

« NEUER DESTABILIZE LARGE FLEXIBLE SPACE STRUCTURES 
WHEN THE ACTUATORS ARE CO- LOCATED WITH THE SENSORS 

e WITHIN THE REDUCED-ORDER DESIGN MODEL, ANY rtlOUNT 
OF DAMPING DESIRED CAN BE ADDED TO ANV MODE EXACTLY 
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NUMERICAL ANALVSI S OF UIBRATION MODES 
t. LOS ERROR DUE TO UNIT INITIAL MODAL DISPLACEMENT 
MODE 1234567 89 10 

PEAK .37 .53 .54 .93 1.3 .14 .51 .002 .18 .03 
= = = = > 5, 4, 3, 2, 7, 1 , 9, 6, 10, 8 


2. MODAL DISPLACEMENT DUE TO RAPID POINTING SLEW 
MODE 123 4 5 

PEAK 21.6 603 41.2 13.7 0.49 

^ 3 * 1, 4, 5, 6, 7, ... 


WHICH MODES REALLY REQUIRE ACTIUE CONTROL? 

NEED AN ALTERNATIUE AND MORE INDICATIUE MEASURE 


3. LOS ERROR SOLELY DUE TO EACH MODE 
BY THE SLEW 

MODE 12 3 4 

PEAK 3.26 88. 6(?) 9.57 6.53 

====> 2 » 3 , 1, 4 (OR 4, 1), 5, 7, 


EXCITED 

5 

0 . 33 
6 , ... 


A 


SOUND MEASURE OF THE SI 

INPUT (SLEW EXCITATION) 
DULY COMBINED 


GNIFICANCE OF EACH MODE: 
AND OUTPUT (LOS ERROR) 
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8.86E+B1 



A 


i I 
t I 


/ 


3.55E-05 , ^v 


LOS ERROR (DEG) DUE SOLELY TO MODE 5 RUN TlflE- 4.892 P f0000.6na 


LCOPY 
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Actuator/Sensor Influence on first 10 Modes 
Mode Act. 1-3 Mode Act. A - 6 Mode Act. 7-8 Mode Act, 9-12 


2 0.30019961E-02 5 0.36487188E-01 

4 0,412203080-03 4 0.25172627E-01 

1 0.401463210-03 3 0.15999462E-01 

3 0.191843690-03 2 0.1559580CE-01 

5 0,111864740-03 7 0.14711439E-01 

6 0.698817890-04 1 0.13037169E-01 

7 0.368294570-04 9 0.57048416E-02 

8 0.262615320-04 6 0.34139471E-02 

9 0.150721070-04 8 0.12352261E-02 

10 0.134977470-04 10 0.63637015E-03 

Mode Sen. 1 - 3 Mode Sen. 4-6 


2 0.143112180+01 2 0.107114020+00 
1 0.140613840+01 1 0.103388680+00 

3 0.81986851E+00 3 0. 10171293E+00 

4 0.39743480E+00 4 0.69439910E-01 

7 0.303959760+00 9 0.68373762E-01 

9 0.25503686E+00 8 0.67025743E-01 
6 0.218527420+00 5 0.63191518E-01 

8 0.146238790+00 10 0.461036000-01 
10 0.108015390+00 6 0.399357790-01 

5 0.743995900-01 7 0.322639120-01 

Mode Sen. 7-8 Mode Sen. 9-12 


2 0.286900670-03 5 0.34890966E-02 2 0.134668560+00 2 0.107114020+00 

4 0.393901280-04 3 0.32387748E-02 1 0.127369450+00 1 0. 10338868E+00 

1 0.391135980-04 4 0.24055073E-02 3 0.124078520+00 3 0.101712930+00 

3 0.189407890-04 2 0.15346858E-02 4 0.381589010-01 4 0.694399100-01 

5 0.106898480-04 7 0.14879148E-Q2 7 0.367935930-01 9 0.683737620-01 

6 0.667790440-05 1 0.15531352E-U2 9 0.308794600-01 8 0.670257430-01 

7 0.351946910-05 9 0.67911280E-03 6 0.208320750-01 5 0.631915180-01 

8 0.250955216-05 6 0.326240990-03 8 0.140055360-01 10 0.461036000-01 

9 0.1/44029380-05 8 0. 118040010-03 10 0,102999060-01 6 0.399357790-01 

10 0.128984480-05 10 0.60812166E-04 5 0.907372120-02 7 0.322639120-01 
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MODAL - DASH POT ItP. 1 


PART 1: LINEAR UELOCITY FEEDBACK GAIN G LUR 

2 FORCE ACTUATOR ON REFLECTORS 
--> U 7 (X AXIS); U 8 CY AXIS) 

2 LINEAR UELOCITY SENSORS AT REFLECTOR END 
"-> V 15 (X AXIS); Y 16 (Y AXIS) 

2 '’MODELED MODES'’ FOR DAMPING AUGMENTATION 
5 1 = 2 X 60 x u»ji = 2.0964 
> TIME CONSTANT^ 0.95 SEC 
MODE 2 : Sj = 2X 67x u»2 = 2.6389 

--> TIME CONSTANT^ 0.76 SEC 


PART 2: ANGULAR UELOCITY FEEDBACK GAIN G* UR 

3 TORQUE ACTUATORS ON REFLECTOR 
--> U 4 CX AXIS); U 5 (Y AXIS); U 6 CZ AXIS) 

3 ANGULAR UELOCITY SENSORS AT REFLECTOR END 

> Yl0 CX AXIS); Yu CY AXIS); Y CZ AXIS) 
3 "MODELED MODES" FOR DAMPING AUGMENTATION 
MODE 3 : 5^=2X3xw 3 = 0.3065 

--> TIME CONSTANT= 6.53 SEC 
MODE 4 : S J = 2X 3x = 0.4470 

MOnp c:. T IME CONSTANT = 4.47 SEC 

MODE 5 . 5 S _ 2X 3x u> 5 = 0.7742 

TIME CONSTANT= 2.58 SEC 
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_ dx 

y = CD Jt 


DYNAMICS : " £? + D |f + „ x = f 

FORCE CTOROUE) actuators and displacement sensors: 

f - Bp u 

CONTROL LAW FOR D I SPLACEMENT- OUTPUT FEEDBACK ; 

u = ~ G D y 

FULL-ORDER CLOSED-LOOP SYSTEM EQUATION: 

d2 n . _ 

dt * A dt + (ft2+ 5 Bp G C D 5) n = 0 

tlOP^ki SPRING ftPPR^ru 

design to augment stiffness to EACH MODE 

OF A REDUCED-ORDER MODEL 
LET “iNEW BE DESIRED FREQUENCY FOR MODELED MODE i 

*et s r Bp g c v 5f1 = diag[„2 new -„2j 

= DIAGjV.l 

then SOLUE FOR FEEDRACK GAIN MATRIX G. 

G = <5 m b f> + DI AGjV.J (C D S M )t 
USING THE PSEUDO- I NUERSES < >t DEFINED AS FOLLOWS 


<5^ Bp 

= <5^ Bp,T j- (5 T Bf) (5 t Bp jT1_, 

<C D s M> f 

= [< C D *M> T <c d <c d b „> 


r -c4 
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PRE-DESIGN ANALVSK.S 


— MODAL SPRINGS 

PLACEMENT OF 2- AXIS PROOF-MASS ACTUATORS 

1 AT REFLECTOR END 

-- PEAK OF MODES 1, 2, 3 

1 AT 92FT FROM SHUTTLE C 70 . 77x LENGTH) 

-- PEAK OF MODE 4 

FFT ANALYSIS OF RPR SLEW DISTURRANCE 

==> SHF IT MODES 2 & 3 UP AND AWAY tf? 

AUOID CONTROL SPILLOUER TO MODE 1 f 
IGNORE MODE 4 
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modal-spring hs.i 

LINEAR DISPLACEMENT FEEDBACK GAIN G LDM 

2 2- AXIS PROOF-MASS ACTUATORS ON MAST: 

1 AT REFLECTOR END 

--> Ug (X AXIS); U 10 (Y AXIS) 

1 AT 92 FT FROM SHUTTLE C70.77* LENGTH) 
~“ >u ll <X AXIS); U 12 CY AXIS) 

4 LINEAR DISPLACEMENT SENSORS ON MAST: 
CO-LOCATED WITH PROOF-MASS ACTUATORS 
~ _> v 13, VjL7 <X AXIS); 
v 14i V 18 (Y AXIS) 

3 ,ri0DE LED MODES* ' FOR STIFFNESS AUGMENTATION 
MODE 1 : (r J = 0 

MODE 2 : erg = C2n X 0.7)2 -(Zn x 0.3136)2 

= 15.4627 

— ° PE 3 : 5 3 = (2n x 0-85)2 - ( Zn X 0.812)2 

= 2 . 4290 
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MODAL-DASHPOT HP. 2 


LINEAR VELOCITY FEEDBACK GAIN G lum 

2 2- AXIS PROOF-MASS ACTUATORS ON MAST: 

1 AT REFLECTOR END 

--> Ug (X AXIS); Uj^ (Y AXIS) 

1 AT 92 FT FROM SHUTTLE C 70 . 77x LENGTH) 
~ _> U H (x AX IS); U 12 <Y AXIS) 

4 LINEAR VELOCITY SENSORS ON MAST: 

CO-LOCATED UITH PROOF-MASS ACTUATORS 
""> Y 15, V 19 CX AXIS); 
y 16* y 20 <Y AXIS) 

3 MODELED MODES’' FOR DAMPING AUGMENTATION 
— 9 PE — = 2X2.7xuii = 0-0943 

MODE 2 : = 2X 2 . 7x w 2 (NEW) = 0.2375 

— PE 3 : 5 3 = 2 x 2.7 *103 = 0.2758 


MODAL-DASHPOT HP. 3 

PART 1: LINEAR VELOCITY FEEDBACK GAIN G LUM 

PART 2: ANGULAR VELOCITY FEEDBACK GAIN G AUR 
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8.82E+0 


2.29E-03 


v A A \ 

V^/v 

LOS ERROR (DEG) DUE SOLELY TO HODE 3 RUM TIHE= 5.000 P f0010.6na 


6.37E+00 


3.45E-03 


!\ i 


A/V \j w 


A A A 

\A a a a / 

a \ . y 1/ v / \/ \/ 


\/V\/Va 


LOS ERROR (DEG) DUE SOLELY TO ItODE 5 RUN TIDE- 5.080 P f0818.6iui 



ETfl3 RUN TIflE = 5.000 P f0010.6na 







6.34E+01 



LOS ERROR (DEG) DUE SOLELY TO MODE 3 RUN TIflE- 4.892 P f0100.6na 
OPY 


05E+01 




2.76E+02 


-2.70E+02 

LCOPY 

L48E+01 


-1.61E+01 




ET05 RUN TINE = 4.892 P f0100.6na 


LCOPY 




5.6BE+01 



f rnnu 






5.48E+01 



7.51E+00 


3.55E-05 


HOP? 




/\/ j\ i\ i 

v I W W V |i 

LOS ERROR (DEG) DUE SOLELY TO flODE 5 RUN TIDE: 1892 P f3100.6 na 
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J.oot+Hl 


LOS ERROR (DEG) DUE SOLELY TO flODE j run T{fJE; 


5.000 P f211B.6na 


! /! 
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I /' i 
. / / 




I fi | pwij/vwu/ 


A U ,\ A A 
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3.39E+01 


-4.36E+00 



VM RUM TME = 5.000 P f3110.6na 
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CONCLUSIONS 


GENERAL: 


® MODAL- DASHPOT AND NODAL-SPRING CONTROLLERS 

PROUIDE QUICK AND EFFECTIUE UIBRATION CONTnm 

•' EUE " excited by most uiolInt; baSg-SaSg 

* Hl S2ZnSI N c£ ROBLEMS CAN BE OUOIDED by 

PROPER ?pu5? T I2 N ° F ,>M0I >ELED MODES'* AND 
PROPER LEUEL OF AUGMENTATION 

R MODAL DASHPOTS AND MODAL SPRINGS MOST EFFECTIUE 

__ M pSn INITIAL PERIOD OF LARGE UIBRATIONS 

PRECIS! oN^pm PERFORMANCE CONTROLLERS FOR 

PRECISION POI NT I NG^STABI LIZATI ON LATER 


9 LOS ERROR DUE SOLELY 
DISTURBANCE PROUIDES A 
OF I ND I U I DUAL MODES 
-- CORRECT SELECTION 


TO EACH MODE EXCITED BY THE 
SOUND MEASURE OF IMPORTANCE 

OF MODES TO CONTROL 


30 ? 


SPECIFIC ON THE HU1ERICAL SIMULATIONS: 


MAST BENDING 


FORCES AND MOMENTS 


* <iRE " T, ' v 

ANn 2^2 UIRE LA *S E CONTROL 

-- 2m-? 221 BE UERY precise 

BUT ARE FAST AND EFFECTIVE 

EXCESS IUE LOS 1 JrrTER R AND MA9 T T r TI ° N PRESENTED 
CF0100 US F0000T ST BENDIMG 

’ us lj;Fr~I SsrsMK'sss ass. 

' 8, ""“ 

effectively and 

(F0110, F2110 


& 


BENDING 
OU1 CKLY 

F3110 US F0000 ) 


SUPPRESSED 


« MORE 


tF3110 US F2110) 

SUPPRESs'f^ 0 ^ C0MTR01 - FORCES 
SUPPRESS LESS LOS JITTER, LESS 


AND MOMENTS, 
MAST BENDING 
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ISSUES NEEDED TO BE ADDRESSED: 

* COUPLING OF RIGID-BODY DYNAMICS 

® INTEGRATED DESIGN WITH LQG/LTR POP ur^u 

-- MODAL DASHPOTS AND SPRINf^ i PRECISION 

72r5^S NCE ST «BVMTY P ^2 G R 0 2gsTf!2Is L °° P 

AS OUTER LOOP TO ENHANCE PRECISION 

8 TO TN\5s M pofN?iJ? E Jn Q ^ I,,I:,> ^CURACY 
ui> POINTING AND STABILIZATION 

» EULUATION ON THE LABORATORY APPARATUS 
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MODAL CONTROL (CONT'D) 

Truncated modal equations: n .(t) + o>?n.(t) = f.(t) + v.(t). 
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MODAL CONTROL (CONT'O) 
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INDEPENDENT MODAL-SPACE CONTROL HMSC) 

Linear (proportional and rate feedback) control: 


i) 


4 
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INDEPENDENT MODAL-SPACE CONTROL (IMSC) (CONT'D 
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REGULATION of 


the scole configuration 


INVESTIGATORS 

Gregory A. Norris 
Emmanuel G. Collins 
Robert E. Skelton 


} 


PERFORMANCE REQUIREMENTS 

(I) MAINTAIN RMS OF THE STEADY STATE LINE-OF-SIGHT (LOS) 
ERROR WITHIN A SPECIFIED BOUND. 

(II) MAINTAIN STEADY STATE ACTUATOR VARIANCES AS CLOSE AS 
POSSIBLE TO SPECIFIED BOUNDS. 

ORIGINAL SCOLE CONFIGURATION 

• LOCATION OF 2. PROOF MASS ACTUATORS NOT SPECIFIED. 

• 42 SENSORS PROVIDED. 
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objectives 


(I) DETERMINE LOCATIONS FOR PROOF MASS ACTUATORS. 

(II) DETERMINE A REDUCED SET OF SENSORS. 

(III) DESIGN A CONTROL LAW TO MEET PERFORMANCE REQUIREMENTS 
FOR LOS ERROR AND ACTUATORS, 

• SOLUTIONS TO THE 3 PROBLEMS ARE INTERDEPENDENT. 

• CHOICE OF ACTUATORS AND SENSORS INFLUENCES CONTROL LAW. 

• CHOICE OF CONTROL LAW INFLUENCES SENSOR AND ACTUATOR 
SELECTION. 
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LINEARIZED DYNAMICAL MODEL 


VECTOR SECOND ORDER MODAL FORM 

• • • 2 _ , 

n + Dn + Jr n * B(u+w) 

output vector y 


y l = L0S x ’ y 2 “ L0S y ’ y 3 = L0S z 
o 2 2 2 1/2 

E(LOS error) = (Ey^ + Ey^ + Ey^) 



measurement vector z 


P,r 


position & rate measurement vector 


where 


z = 


P n + p n + v 
P V p,r 

acceleration measurement vector 

• • 

Qn + v 

a 

2 • _ _ 

“* Dr) + Bu + Bw) + v 
' a 


z = 


P*r 

z -QBu 


M n + M i) + v 
P v 


M = 
P 


M 


0 


L 0 -qjt| 

- p o' 

v 

0 -QD 
v 


P*r 


. v +QBw] 

La 1 


J 

ASSOCIATED SENSOR NOISE (v) & ACTUATOR NOISE (w) ARE 
CORRELATED 
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• MODEL OBTAINED USING CUBIC BEAM ELEMENT SHAPE FUNCTIONS 

FOR BEAM BENDING AND LINEAR SHAPE FUNCTION FOR BEAM 
TWIST. 

• 32 MODES IN ORIGINAL MODEL. 

• MODAL COST ANALYSIS USED TO REDUCE TO 23 MODE DESIGN AND 
EVALUATION MODEL. 
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modal mode 
cost rank no. 


MODAL COST ANALYSIS 

modal cost fi eq. <hz) 


mod € 
type 


1 

1 

infinite 

0 

rigid body 

2 

2 

infinite 

0 

rigid body 

3 

3 

infinite 

0 

rigid body 

4 

5 

, 91 1 c-+07 

. 299e+00 

bending < rol 1 ) 

5 

7 

. 363o+07 

. 1 18ei01 

bending 

6 

4 

, 336e+07 

. 2760+00 

bending (pitch) 

7 

6 

. 1380+07 

. 81 le+OO 

bending 

8 

8 

, 955e+06 

. 205e + 01 

bending 

9 

10 

. 673C-+04 

. 551 e+01 

b c n d i n y 

10 

9 

. 556c i 04 

. 478ch 01 

bending 

11 

11 

. 246e*i 02 

. 123e+02 

bending 

12 

14 

. 365e+01 

. 24 3 e +02 

bending 

13 

17 

. 24 5e+C 1 

. 393e+02 

twi ?t 

1 A 

12 

. 305e+0Q 

. 129e * 02 

bending 

1 5 

1 6 

. 1 1 6e +00 

. 390e+02 

b en ding 

16 

if. 

. 349e-01 

. 256e+Q2 

bend) no 

17 

26 

. 99 5 e- 02 

. 1 09e + O3 

'J 

b en d i n y 

18 

25 

. 377o-02 

. lC'3e + 03 

b c Ti cl i n u 

19 

13 

, 3760-02 

, 237o h 02 

bend ir, o 

20 

29 

. 174o-02 

. 140ei03 

b i n d i n u 

PI 

35 

. B 36c -03 

. 2 1 5e + 03 

i. e ri d i T i q 

PP 

20 

. 597o-03 

. 586e+02 

b L- Vi c i n g 

23 

28 

. 370 e- 03 

. 1 39o-* 03 

bending 

PA 

23 

. i 25e-03 

. Lil7e-i 02 

bend i ti g 

25 

19 

. 3 1 Oo -04 

. 5810102 

bending 

26 

34 

. 2750-04 

. 2 1 5-t i 0 J 

bend i n q 

27 

32 

. 6170-05 

. 173gi03 

bending 

28 

31 

. 2940-05 

. 175e+03 

bending 

29 

27 

. 131 e-05 

. 135C+03 

twist 

30 

24 

. 140e-07 

. 106ei03 

t li'iGt 

31 

30 

. 1 340-07 

. 1 67 f +03 

tu-i st 

32 

33 

. 4130-03 

. 2C>0t-i 03 

t u» i r. t 

33 

22 

. 2980-10 

. 81 I 01 O 2 

bending 

34 

18 

. 3400-1 1 

. 51 5ci 02 

tu'i ?. t 

35 

21 

. 226o-13 

. 782 ei 02 

twist 

• 

FIRST 5 

FLEXIBLE MODES 

DOMINATE 

MODAL COST 


• BEAM BENDING DOMINATES MODAL COST 
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CONTROL LAW DESIGN VIA 
THE OUTPUT VARIANCE ASSIGNMENT ALGORITHM 

• ITERATIVE ALGORITHM DEVELOPED BY SKELTON AND DELORENZO 

• OBJECTIVE IS TO CHOOSE DIAGONAL Q AND R IN THE LQG COST 
FUNCTIONAL 

v = E^C y T Qy + u T Ru) 

S.T. THE LQG CONTROL LAW SATISFIES 

38 (or < a^) V 1 - 1 n 
WHILE MINIMIZING 


n 

r. 2 

u 

E u. 

E 

1=1 

o» i 

u 2 • 

V. 


bounds on input variances 
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SENSOR AND ACTUATOR SELECTION 
VIA INPUT/OUTPUT COST ANALYSIS 


SUBOPTIMAL APPROACH. 


BASED ON DECOMPOSING COST FUNCTION 


as 


v = E Jy T Qy + u T Ru) 


n 


n 


u 


u 


v = Z v' + z v 
i=l 1 i=l 1 


n 


n 


u z 

v - E v” + Z v v . 
i=l 1 1=1 1 

0 DEFINES ACTUATOR EFFECTIVENESS, 


»r‘ ■ »; - 

AND SENSOR EFFECTIVENESS 


v sen « v v 
i i * 


DELETES ACTUATOR(S) OR SENSOR(S) WITH LOWEST 


EFFECTIVENESS VALUES. 
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SOLUTION PROCEDURE 

BEGIN WITH LARGE SET OF PROOF MASS ACTUATORS 
FIXED LOCATIONS 


► 

DESIGN LQG CONTROLLER 


1 


delete actuators 

Ki Q 

i 

*1 w 

final no. of actuators? 

1 

1 yes 

— i 

DESIGN LQG CONTROLLER 


j 

i 

DELETE SENSORS 

mo r 

HZn~ 

n r 4 ■* * « 1 


1 FINAL NO. HP d?KTC*An r* r% 1 

I w w* uunovivo : j 

Tves 


DESIGN FINAL CONTROLLER 








SOME RESULTS 


ORIGINAL SCOLE PROPOSAL 

rms(los error) < .02 deg 

OUR findings 

if noise through shuttle cmgs only: 
nns(los error) > .045 deg 
if equivalent noise through all actuators: 
nns(los error) > .075 deg 

CONCLUSIONS 

• ORIGINAL SPECS ON LOS ERROR ARE NOT ACHIEVABLE. 

• MUST MODIFY LOS SPECS. 
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ACTUATOR SELECTION 


DEFINE 


USUM - 


"u EuJ 

1 2 

1=1 


dimensionless measure of 
total control effort 



FINDINGS 

• BY USING REDUCED SET OF ACTUATORS THERE IS A 50% 
SAVINGS IN CONTROL EFFORT (AS MEASURED BY USUM). 
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• PROOF 

beam 


AQ-Tu^ToR^ NE/VR 

ARE More E ff e- clt i o el 


\ o 9 


of rne 
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"Sensor ^eleclt\oh 



NO. OF 


• GOOD PERFORMANCE MAY BE ACHIEVED WITH A MUCH SMALLER 
SET OF SENSORS, 
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CONCLUSIONS 


(I) achievable R0R) * - 02 DEG IS N0T 

RMS (LOS ERROR) — .05 DEG IS 

ACHIEVABLE IF NOISE IS ONLY THROUGH 
SHUTTLE CMG'S . iHKOUGH 

RMS(LOS ERROR) Si. .08 DEG IS 
ACHIEVABLE IF (EQUIVALENT) NOISE 
IS THROUGH ALL ACTUATORS. SE 

(II) p ??™® 35 ACTUATORS SHOULD BE 
PLACED NEAR TOP OF MAST. 

(Ill) GOOD PERFORMANCE MAY BE ACHIFVPn 

WITH A (SIGNIFICANTLY) REDUCED^ 

SET OF SENSORS. REDUCED 
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Evaluation of On-Line 
Pulse Control for 
Vibration Suppression 
in Flexible Spacecraft 

by 

G. A. Bekey 

S. F. Masri 

R. K. Miller 

Univ. of So. California 
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• h 

EVALUATION 

PULSE CONTROL FOR VIBRATION 
SUPPRESSION IN FLEXIBLE 
SPACECRAFT 

G.A.Beketf. S.F. Masri R. K. Miller 

University of Southern California- 

Los Aneeles, C\ 

OUTLINE 
l. Objective 

E. Modeling Issues 
—Beam vs. Truss 
-NL- FEM t numerical problems 
M. Control Issues 

—ED TLdlse \ctuator Development 
-Pseudo "Pulse AJteoritAm Dev. 
—Large NL Simulation "Problems 
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OF ON-LINE 




OVERALL OBJECTIVE 



0-00970 

roA/s 2 - 


-0 009*0 
md./s 1 ’ 



(.Gent /Un, (984 SCDLE Mia. ) 



Minimum Time 20 * 

Stew Maneuver 


Mess-Ejection Rise Gontret Strategy. : 


P,(t) 


-c i »9n(v,)|vj 


0; 


V < ' < <W 

i i i 

(t n ♦T. ) < t < t n 
°, d i Vi 


Typical Experimental Results : 
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MODELLING ISSUES 

l Continuous Beam vs. Truss 

* AtieX / TorsioneJL Coupling 

* Local Member fhrtCcuoo tion in Mo 

* Parametric Resonance Roblams 





Mast FMoAt Beam 
Alternating Bay Identical Bay 

Diagonals Diagonals 
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LINEAR TRUSS RESULTS 


Z. Linear finite Element Model Characteristics 
•COFS-l Hardware Configuration 
•54 Bays, 60m 

• 171 nodes, 466 elements, 522 V.of F. 

•July 1986 data, for member characteristic 5 
from Aftro Aerospace Carp. /Harris Corp. 
•Match modal results with Astro /Harris 
•Transient Response Simulations : 
-Rayleigh damping : f 7. t 107» 

-Sine -sweep, tip exctation 
— NonStetionaru Random, tip excitation. 

— Harmonic , base excitation 
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FIG. 1. NONLINEAR THREE-DIMENSIONAL FINITE ELEMENT MODEL OF COFS I MAST 
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34 - 4 


FIG. 4. REPRESENTATIVE MODE SHAPES CORRESPONDING TO THE LINEARIZED VERSION OF THE USC 
NONLINEAR COFS I MAST FINITE ELEMENT MODEL. 




ACCELERATION EXCITATION 



FILE.-NAC9N171 



TIME SEC 


Sivept Sine "Response - Tip Excit 

FILE:NAC9N171 
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Nongfctionaru Random Response- Tip Excii 







NONLINEAR UNITS ELEMENT MODEL 


Type of Nonuntarity 


-e 


(f )* 2 * 

0-57\f\[$-int,.t4f)1 

0£ If I <0-63 

*0-86; If I >0-63 
Hert&ian Joint Contact 
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NONLINEAR SIMULATION PROBLEMS 

1- Excessive CPU Time 

• 92+ Hours ior T* 3 Fundamental Periods 

on VAX H/750 

•Small At Tjj / looo) required for 

numerical Stabilitu 

2. Modal Order Reduction Necesseru 
•Non parametric "RONN" Modal - in process 

• farnmetric /Super element Model- in process 
• Validity ??? 



352 



CONTROL ISSUES 


l Simulation 


PivoresS 



* 

■inppi 

•Lra. Displ. 


*NL, Adaptive, 

•Joint Nontin. 


’Intermittent 

• R- K, NewmarkB, 


’ Special Case: 

Wilson Q t etc . 

1 

•i 

,? 


1-Tune step interupt redd. • 

2 . Excess Storage for State variables 

3 . Stability and. restart, capability 
jar tune- Pepping algorithms 
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THREE-DIMENSIONAL FINITE ELEMENT MODEL OF COFS I MAST 
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ACTIVE STABILITY AUGMENTATION OF LARGE SPACE STRUCTURES: 
A STOCHASTIC CONTROL PROBLEM 


A.V. Balakrishnan 

Department of Electrical Engineering 
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Paper to be presented at the IFAC Conference 
Vilnius, May 1986. 


on Stochastic Control, 


Research supported in part under NASA grant no. NAG 1-464. 
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1/ Introduction 


In the 1987-1990 period NASA has planned several ground and flight 
experiments with the eventual objective of deploying large flexible structures 
in space. A currently active precursor is the SCOLE experiment [1]. Here 
the problem is that of slewing an offset antenna on a long (130 ft.) flexible 
beam-like truss attached to the space shuttle, with rather stringent pointing 
accuracy requirements (±.02 degrees). This paper examines the relevant 
methodology aspects in robust feedback-control design for stability augmen- 
tation of the beam using on-board sensors. We frame it as a stochastic 
control problem - boundary control of a distributed parameter system described 
by partial differential equations. While the framework is mathematical, the 
emphasis is still on an engineering solution. 

The fact that the deployment is in space makes model uncertainty the 
major consideration in control design. Particularly serious in this regard 
is for instance the modelling of inherent damping in the system long known 
to be difficult [2] , and a still unresolved problem even in theory. Hence 
robustness becomes a must feature, even at the expense of optimality. Another 
aspect is the complexity of computation, making any simulation study a costly 
undertaking. 

The overall model involving both slewing and beam stabilization is still 
not well understood. Hence the two problems — of slewing and stabilization — 
are best studied, at least in initial efforts such as reported here, separately. 
We attempt stabilization at the termination of the slewing so that in parti- 
cular the system is essentially linear except for a small nonlinear term 
contributed by the kinematic nonlinearity. It should be noted that at present 
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we do not have a stochastic time-optimal control theory adequate for optimal 
slewing based on sensor data. 

An abstract mathematical formulation is developed in Section 2 as a 
nonlinear wave equation in a Hilbert space. We show that the system is 
controllable and develop a feedback control law that is robust in the sense 
that it does not require quantitative knowledge of system parameters. The 
stochastic control problem that arises in instrumenting this law using 
a PP ro P r iate sensors is treated in Section 3. Using an Engineering first 
approximation which is valid for "small" damping, formulas for optimal choice 
of the control gain are developed. 

2. Abstract Formulation 

We are concerned with the mast stabilization problem only and the model 
we use assumes that the angular velocity of the shuttle-antenna system is 
small enough to be neglected. We model the mast as a thin prismatic beam. 
There is then the question of whether a finite-element model or a continuum 
(involving partial differential equations) model should be used. Here we 
deal only with the latter, the basic governing equations being beam bending 
and torsion equations with controls at the boundaries. 

With reference to Figure 1, the beam of length L is along the Z 
axis, z being zero at the shuttle end. u^(-), u 0 (*) will denote the 

displacements along the Y-Z, x-z planes and u^(*) the angular deflection 
about the Z axis. In addition proof-mass controllers are provided at 
points s 1 and s 2 , on the beam, the locations to be chosen optimally. 
Control moments are applied at both ends as well as control forces at the 
reflector center. The various moments of inertia and masses are specified 
in [1], (2]. 
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E 

FIGURE 1 

SHUTTLE/ANTENNA CONFIGURATION 
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We first develop an abstract mathematical model. We define 
H - L 2 rO,L ) 3 x R 14 o < L < co 

with the usual inner-product thereon denoted [ , ]. We fix the points 

° < s 2 < s 3 < L and define a linear operator A into H with domain D 
m H defined as follows. We use u^(-), u 0 (-), u^-) to denote the 

functions in L_,[0,L] 3 . Thus an element x in H is denoted 

V ,} 

v*> 




X 17 



The domain 

D consists of elements 

x such 

that u', 
<P 

u e' € V°' L] 

and u'" ( • ) 

has - derivatives in 

[0,s 2 ] , 

fS 2 ,S 3 1 

and [s 3 ,L]; 


similarly for u Q ( • ) ; ( • ) such that u^(*) and uj(-) « L 2 [0,L] ; the 

remaining components of x are specified as 


X 4 

“ 

u. (0+) 

V 

X 11 

= 

u;(L-) 

V 

X 5 

= 

u e (0+) 

X 12 

= 

u e (L_) 

X 6 

= 

u.(l-) 

X 13 

s 

V L_) 

X 7 

= 

u e ( L-) 

X 14 

s 

% (s 2 ) 

X 8 

= 

u! (0+) 
4 1 

x is. 

s 

u e (s 2 ) 

X 9 

= 

u 0 (0+) 

x l6 

- 

W 

X 10 

s 

V 0+ > 

x,. 

m 

u»(s.) 


Thus at least for x in D, we may identify the finite-dimensional part 
as the "boundary." The operator A is then defined by 

y * Ax 

where the functional part (in L 2 (0,L] 3 ) is given by 

ei u;-(.) 

9 ♦ 

EI e u e" ( * ) 

-GI^uJC) 

and the boundary part by: 


y 4 

X 

El .u"' (0+) 
9 ♦ 

y n 

x 

EI.u"(L-) 
9 9 

y 5 

s 

El gUg" (0+) 

y l2 

- 

E W L-) 

y 6 

*= 

-El u *" (L-) 
♦ $ 

y l3 

X 

G V* (L " ) 

y 7 

= 

- EI e u e* (L -> 

y i4 

- 

E V u i ,,(s 2 +) - u ; ,,(s 2- n 

GO 

>1 

- 

-El u"(0+j 
9 9 

y i5 

= 

EI e (u e' (s 2 +) - u 9" (s 2‘ )) 

y 9 

- 

" EI e u e (0+> 

y i6 

* , 

E V u r ( v> - u ; ,,(s 3’ )) 

y io 

* 

- G V* (0+) 

y !7 

m 

E1 e (u e” (s 3 +) ‘ u e ,(s 3 - )) 


It may then be verified that D is dense and A is self-adjoint and 
nonnegative definite. Moreover A has a compact resolvent with a complete 
orthonormal set of eigenfunctions (modes). Zero is an eigenvalue. 

The control system dynamics can then be characterized as a nonlinear 
wave-equation: 

Mx(t) + Ax ( t) + K (x ( t) ) + Bu (t) = 0 (2.1) 

where M is a 17 * 17 nonsingular nonnegative definite matrix, and defines 
self-adjoint positive definite linear operator H onto H. The control 

u(*) is in R 12 , and 
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Bu 


X - co1 - (0 ' U 1' u 2 , u 3> u 4 , u 5 , u 6 , u 7 , u 8 , u 9 , u 10< u xl , u 12 I 

We have thus only "boundary" control. The nonlinearity is kinematic: 


where 


K(x) 


0 

0 

0 

0 

0 

0 

0 


n ! • 'A 

n 4 • V. 


0 

0 

0 

0 


n l = 

col 

(x e ,x 9 ,x 

ft 

4 

col 

(X 11' X 12 


^4 are symmetric positive definite (moment) matrices and 
• denotes vector cross-product. 

Two relevant properties of the function K(-) are: 

(i) [ K (x) , x) = 0 

(ii) 1 K ( x ) || < const. ||x|| 2 

We do allow for "state noise" and let 


4 


A 
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A , and the 


N(t) 


FN (t) 


N x (t) 
N 2 (t) 
N 3 (t) 


x ( t) 


where N(t) is white Gaussian with spectral density matrix 
components of x(t) are defined by 


i = 1 f • • • f 7 


x i (t) = 0 

Xg(t) = ^(t) 

x 9 (t > = N 2 (t) 

X 10 (t) " N 3 (t) 

x.(t) =0 i > 10 . 

Note that the "boundary" values are part of the state. 

State-space Form 
With 


Y(t) = 


x(t) 
x (t) 


we go over to the state-space form: 

Y(t) = Ay( t ) + K(Y(t)) + 8u(t) + F(N(t)) 

where 


8u(t) 


-M -1 A 0 


-M _1 Bu(t) 


( 2 . 2 ) 
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and in the notation 


we have 


Y = 



Y c // x // 


ICC Y) 


FN(t) 


0 

-M -1 K(y 2 ) 

0 

-M -1 FN (t) 


well known, we can introduce a new inner product, the "energy" 
inner product 

[Y,Z] £ = [7a y x , ✓a zj + [My 2 , z,,) 

on /?(A) x « . i?(A) is the orthogonal complement of the null space of A. 

We denote the completed space by We shall from now on consider only 

We have: 

A + A* = o 

and of course A has a compact resolvent and we have an orthogonal decompo 
sition of H given by 


Y = 



where P R is a two-dimensional projection for each k, 



A Vk 


( 2 . 3 ) 


spanned by 
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where 




*K • 


2 ^ _ 2 

\ > 0* <*> k •; (2. 


tM* k , ♦ J 


‘ 5 - 


Let S(t) denote the semigroup generated by A. Then we have the 
sentation: 


repre- 


S (t ) Y = 


l S(t) P Y 
1 K 


More explicitly, if 


Then 


P k S(t) P k 


S (t) Y = 


s( t) P R . 


y x (t) 


y 2 (t) 


y 2 (t > = y 2 (t) 


and 


CO 

V l (t) * l fy l' M *k ] *k COS V + [ ty 2' M V * k oT 


sin 


V 


Note that it is required that y satisfy: 


(2.5) 


l ty r M# k ] 0,2 < 


It is easy to establish existence and uniqueness of solution for the 
integral version of (2.2): 

t t 

Y ( t) * S (t) Y (0) + / S(t-o) 8u(o) do + / S(t-o) Fn(o) do 

0 0 

t 

+ / S(t-o) K(y (o ) ) do , (2.6) 

0 

without invoking any nonlinear semigroup theory, by just Picard iteration. 
See [3] . 
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We can now state the basic result that yields a robust feedback-control 
law for the deterministic system (seeing F = 0) . 


Theorem 2.1. 


Let P be any 12 * 12 symmetric nonnegative definite nonsingular matrix. 
Then the feedback control 


u(t) = -P 8* Y( t ) 

is such that the "closed-loop" system 

Y(t) = Ay(t) - 8pB*Y(t) + K(Y(t)) 
is globally asymptotically stable. That is to say 


l*<t)fl E + 0 


(2.7) 


( 2 . 8 ) 


as t 00 

Proof. We refer to [4] for a proof. The proof exploits the fact that (A, 8) 
is controllable in an essential way. in particular the 
generated by (A - 8 P 8*) is strongly stable: that is to say: 


semigroup S (t) 

B 


|S B (t)Yfl E ■* 0 


as t + 00 . 


We also obtain that 


l ( p 8 *Vt)Y, 8 *s B (t)Y) dt = | ||y||2 . 

The control law is also optimal for the quadratic cost functional: 

°° OO 

/ I|/F B*Y(t) II 2 dt + / llu(t) II 2 dt 
0 0 


for the linear system 


Y(t) = Ay (t) + 8/p u (t) 


(2.9) 
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3. Stochastic Control 


To instrument the control law 

u (t) - PB*x(t) 

- Pb(t) . (3.1) 

We need to assume co-located (rate) sensors. The sensor output v(t) 
would then be: 

v(t) = b(t) + N Q (t) (3.2) 

where N Q (t) represents the sensor noise, modelled as white Gaussian with 
(12 * 12) spectral density matrix D. In terms of the state-space represen- 
tation (2.2), we can rewrite (3.2) as 

v(t) = CY ( t) + N Q (t) (3.2) 

where 

C = B* 

and C is of course finite-dimensional. If we assume that the separation 
principle applies, a reasonable choice of control law would be 

u(t) a Pb(t) (3.3) 

where, E denoting conditional expectation: 

b(t) = E[b(t) j v ( s) , s £ t) 

and of course 

b(t) = CY ( t) 

a 

where Y(t) is the Kalman state estimates 

Y(t) = E[Y(t) | v (s) , sit] . 
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Even if we were to neglect the nonlinear term K(-), this would require an 
infinite-dimensional Kalman filter, which even if we could instrument it, 
would depend on quantitative knowledge of the system parameters. Hence this 

filter would need to be simplified in considerable measure, in favor of 
robustness . 

The simplest version would be one that did not distort b(t) and thus 
would lead to the control law: 

u(t) = Pv(t) . (3.4) 

We are thus introducing a noise input into the system which may excite 
higher-order modes. Let us therefore study the system response which is 
now given by the stochastic equation: 

Y(t) = (A - 8p8*) y ( t ) - Bptyt) + K(y ( t) ) + Fn ( t) (3.5) 

This can be expressed as an integral equation: 

t 

Y(t) = Y 0 (t) + / S B (t-o) K( Y(o)) do (3.6) 

where 

t t 

Y o (t) = s B (t ) Y < 0 > - / s (t-o)BPN (o)do + / s (t-o) F n(o) do 

0 u 0 B 

(3.7) 

We note that because K( • ) is locally Lipschitzian , we may solve (3.6) by 
Picard iteration: 

t 

Y n +1 = Y o (t) + / S B (t-o)K(Y n (o)) do . (3. 8) 

We omit the details; see [ 3 ]. More important to us is actually (3.7). 

We want to show that the process Y^-) is asymptotically stationary and 
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evaluate its covariance function. Following [5], since S B (*) is strongly 

stable, it is only necessary to show that for Y in Hi 

E 


/ is (0)Bpdp8*S.(0)*Y, Y]_ do < 

o B B E 


and also that 


/ [S (o)FaFV (o)*Y, YJ _ da < 
0 B B E 


(3.10) 


For this purpose we note that Sg(t)* is strongly stable with generator 

A* - 8p8* 


and analogous to (2.9) we have that 


J ll^e*S n (t)*Y|| 2 dt = i |y|| 2 . 


Hence 


Since 


oo 

/ | Sd p8*s (t)*y|| 2 dt J 


< - . 


I f *s b (t) * y| < |8*s B (t)*Y| 


we also obtain (3.10) . For Y, Z in H let 

E 


[R(t,s)Y, Z] = E ( [Y q (t) , Y ] [ Y q ( s ) , Z ] ) . 


Then we have that 


R(t,s) = S (t-s) R(s,s) , 

D 


t £ s 


and bonce it follows that 


limit R (t+L, s+L) « S (t-s) R , t i s (3.11) 

L-*® a " 


* / |/Drt*S (t)*v| 2 dt ♦ / |/AF*S.(t)*V| 2 dt . (3.12) 
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The process Y Q ( * ) is thus asymptotically stationary with covariance operator 

S B (t ' s)R « t * s 

note that is not necessarily nuclear, even though R(t,t) will be if 

R(OfO) is. Indeed taking 


D = dl ; 


P = I 


we obtain that 


/ | /d pB*s b (t) *y| 2 dt = | JyK 2 . 

0 2 E 


From (3.8) we can show that the process Y(t) is asymptotically stationary, 
since y.) will have this property for each n. Since it would appear 
that the nonlinearity is small, we shall now concentrate our attention on 

the linear approximation Y ( • ) . 

o 

The eigenfunctions of (A - BpB*) * ^ 

) are approximately the same as that 

of A and the eigenvalues are 

a k 


°k 1 iw k ' 


« 1 


* where 


Hence 


2 °x - [pb k' b k ' • 


+ + , 


lR,f f, , J (IDPb k' Pb K^I-*,.r-K]) 
k k E k [Pb k' b k ] 


(3.13) 


(3.14) 


which is thus the noise energy in the kth mode. We see that increasing 

the damping but also increases the noise excitation. In prac- 
tice one would want a compromise between increasing damping at selected 
low order modes but keeping the noise excitation ,t higher order modes within 
hound. Clearly further work is needed before any attempt at control design. 
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We may also mention one point of purely theoretical interest. To 
characterize the distributions of the noise response of a nonlinear system 
described by ordinary differential equations one uses the Fokker-Planck- 
Kolmogorov equations which are partial differential equations. In (3.5) 
we have a nonlinear partial differential equation; it would be of interest 
to develop a corresponding tool to study the distributions. 
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SUMMARY 

The problem of controlling large, flexible space systems has been the 
subject of considerable research. Many approaches to control system 
synthesis have been evaluated using computer simulation. In several cases, 
ground experiments have also been used to validate system performance under 
more realistic conditions. There remains a need, however, to test 
additional control laws for flexible spacecraft and to directly compare 
competing design techniques. In this paper an NASA program is discussed 
which has been initiated to make direct comparisons of control laws for, 
first, a mathematical problem, then an experimental test article is being 
assembled under the cognizance of che Spacecraft Control Branch at the NASA 
Langley Research Center with the advice and counsel of the IEEE Subcom- 
mittee on Urge Space Structures. The physical apparatus will consist of a 
softly supported dynamic model of an antenna attached to the Shuttle by a 
flexible beam. The control objective will include the task of directing 
the line-of-sight of the Shuttle/antenna configuration toward a fixed 
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target, under conditions of noisy data, limited control authority and 

random disturbances. The open competition started in the early part 

of 1984. Interested researchers are provided information intended to 

facilitate the analysis and control synthesis tasks. A workshop is planned 

for early December at the NASA Langley Research Center to discuss and 
compare results. 


INTRODUCTION 


y future spacecraft will be large and consequently quite flexible. 
As the size of antennae is increased, the frequencies of the first flex- 
ible modes will decrease and overlap the pointing system bandwidth. It 
will no longer be possible to use low gain systems with simple notch 
filters to provide the required control performance. Multiple sensors and 
actuators, and sophisticated control laws will be necessary to ensure 

stability, reliability and the pointing accuracy required for large, 
flexible spacecraft. 

Control of such spacecraft has been studied with regard given to 
modeling, order reduction, fault management, stability and dynamic system 
performance. Numerous example applications have been used to demonstrate 
specific approaches to pertinent control problems. Both computer simula- 
tions and laboratory experiment results have been offered as evidence of 
the validity of the approaches to control large, flexible spacecraft. 
Concerns remain, however, because of the chronic difficulties in control- 
ling these lightly damped large-scale systems. Because of these concerns 
and because of the desire to offer a means of comparing technical 
approaches directly, an NASA/IEEE Design Challenge is being offered. An 
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experimenCal test article is being assembled under the cognizance of the 
Spacecraft Control Branch at the NASA Langley Research Center with the 
| advice and counsel of the IEEE (COLSS) Subcommittee on Urge Space 

Structures. This Spacecraft Control Laboratory Experiment (SCOLE) will 
t serve as the focus of a design challenge for the purpose of comparing 

directly different approaches to control synthesis, modeling, order 
reduction, state estimation and system identification. 

The configuration of the SCOLE will represent a large antenna attached 
to the Space Shuttle orbiter by a flexible beam. This configuration was 
chosen because of its similarity to proposed space flight experiments and 
proposed space-based antenna systems. This paper will discuss the "Design 
Challenge" in terms of both a mathematical problem and a physical experi- 
mental apparatus. The SCOLE program is not part of any flight program. 

SYMBOLS 

a acceleration vector ft/sec^ 

A beam cross section area 

c observation matrix 

d noise contaminating direction cosine matrix measurements 

e line-of-sight error 

E modulus of elasticity 

f concentrated force expressions 

F4 force vector 

g concentrated moment expressions 

„ GI torsional rigidity 

I moment of inertia matrix for entire Shuttle/antenna configuration 
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moment of inertia matrix. Shuttle body 

moment of inertia matrix, reflector body 

beam cross section moment of inertia, roll bending 

beam cross section moment of inertia, pitch bending 

beam polar moment of inertia, yaw torsion 

length of the reflector mast, beam 

control moment applied to the Shuttle body 

control moment applied to the reflector body 

disturbance moment applied to the Shuttle body 

mass of entire Shuttle/antenna configuration 

mass of Shuttle body 

mass of reflector body 

mass density of beam 

beam position variable 

direction cosine matrix, Shuttle body () = T () 

earth 1 Shuttle body 

direction cosine matrix, reflector body () = T () 

earth 4 reflector 

inertial velocity. Shuttle body b ° dy 

inertial velocity, reflector body 

lateral deflection of beam bending in y-z plane 

lateral deflection of beam bending in x-z plane 

angular deflection of beam twisting about z axis 

position variables 

displacement of proof-mass actuator 

line-of-sight pointing requirement 

noise contaminating angular velocity measurements 


t 


1 
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0,<t»,y pitch, roll, heading 

C damping ratio 

4 t noise contaminating acceleration measurements 

“l angular velocity of Shuttle body 

t o >4 angular velocity of reflector body 

DISCUSSION 

The objective of the NASA- IEEE Design Challenge concerning the control 
of flexible spacecraft is to promote direct comparison of different 
approaches to control, state estimation and systems identification. The 
design challenge has principal parts, the first using a mathematical model, 
and the second using laboratory experimental apparatus. The specific parts 

of the Spacecraft Control Laboratory Experiment (SCOLE) program will be 
discussed in detail. 

Control Objectives 

The primary control task is to rapidly slew or change the line-of- 
sight of an antenna attached to the space Shuttle orbiter, and to settle or 
damp the structural vibrations to the degree required for precise pointing 
of the antenna. The objective will be to minimize the time required to 
slew and settle, until the antenna line-of-sight remains within the 
angle 6. A secondary control task is to change direction during the 
"on- target" phase to prepare for the next slew maneuver. The objective is 
to change attitude and stabilize as quickly as possible, while keeping the 
line-of-sight error less than 6. 
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Math Model Dynamics 


The initial phase of the design challenge will use a mathematical 
model of the Shuttle orbiter/antenna configuration. It is necessary to 
obtain a balance, of course, between complex formulations which might be 

more accurate and simplified formulations which ease the burden of 
analysis. 

The dynamics are described by a distributed parameter beam equation 
with rigid bodies, each having mass and inertia at either end. One body 
represents Space Shuttle orbiter; the other body is the antenna reflector. 
The equations for the structural dynamics and Shuttle motion are formed by 
adding to the rigid-body equations of motion, beam-bending and torsion 
equations. The boundary conditions at the ends of the beam contain the 
forces and moments of the rigid Shuttle and reflector bodies. The 
nonlinear kinetmatics couples the otherwise uncoupled beam equations. 
Additional terms represent the action of two, 2-axis proof-mass actuators 
at locations on the beam chosen by the designer. 

The rigid-body equations of motion for the Shuttle body are given by: 

“1 ' - ^‘‘vri + ", * % * y,> 
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Similarly, for the reflector body. 


\ ■ ' * "4 + N,4> 


F, + F„ . 
4 B.4 


The direction cosine matrices defining the attitudes of the Shuttle and 
reflector bodies are given by: 


•T ~ T 
T = - u) T 

1 Vl 


*T ~ T 

T , « - w, T, 
4 4 4 


The direction cosine matrices defining the attitudes of the Shuttle and the 
reflector bodies are related to the beam end conditions. 



1 

0 

0 


cosAO 

0 

sinAG 


cosA^F 

-sinAY 

0 

T 4 = 

0 

cosA<)> 

-sinA<j> 


0 

1 

0 


sinAf 

cosAf 

0 


0 

sinA4> 

cosA<£ 


-sinAG 

0 

cosAG 


0 

0 
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where: 
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The equations of motion for the flexible beam-like truss connecting the 
reflector and Shuttle bodies consist of standard beam bending and torsion 
partial differential equations with energy dissapative terms which enable 
damped modes with constant characteristics for fixed, though dynamic, end 
conditions. The system of equations can be viewed as driven by changing 

end conditions and forces applied at the locations of the proof-mass 
actuators. 

ROLL BEAM BENDING: 


a 2 

3 u x 


PA 


3t 


- 2C. /PA El . 


.3 

3s 2 3t 


aV 


+ El 


♦ - n f* (s-s ) 

y 3s n =l ^ ,n n ‘M 3s n 


PITCH BEAM BENDING : 
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3 u,. 


PA 


3t 


r ~ 2 S 
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as^t 0 


aV 


3s 


+ ®0,n H 


YAW BEAM TORSION: 
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a u„, 
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(SHUTTLE BODY FORCE} 
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{PROOP-MASS ACTUATOR} 
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{REFLECTOR BODY FORCE} 
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{SHUTTLE BODY, MOMENTS} 


g 0,2 = 0 


(PROOF-MASS ACTUATOR, MOMENT} 


g 0,3 = 0 


(PROOF-MASS ACTUATOR, MOMENT} 


g s .& 1 " V 4 + VA + M 4 + Vb,^ 


{reflector body, moment} 


The angular velocity of the reflector body is related to the Shuttle 


O 13 c O 
R b = -130 o o 
o o o 
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The line-of-sight error described in figure 2 is affected by both the 
pointing error of the Shuttle body and the misalignment of the reflector 
due to the deflection of the beam supporting the reflector. The line-of- 
sight is defined by a ray from the feed which is reflected at the center of 
the reflector. Its direction in the Shuttle body coordinates is given by: 


where : 


LOS 



[ R I<\ - v • ■>*] 

IK - “p- 2 

r a r r ~ r f 

• r a| 


Rp is the feed location (3.75, 0, 0) 

Rr is the location of the center of the reflector (18.75, -32.5, 

-no) m an undeflected state. 

R A Is a unit vector in the direction of the reflector axis in 
Shuttle body coordinates 


The vector R A can be related to the direction cosine attitude matrices 
for the Shuttle body, Tj , and the reflector body, , by 



The relative alignment of the reflector to the Shuttle body is given by 
T 

T 1 T 4 which is a function of the structural deformations of the beam. 
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The line-of-sight error, e, is the angular difference between the 
target direction, given by the unit vector, D T , and the line-of-sight 
direction in Earth axes, T^qs* 

ARCSIN |D T X T iRl()S | or ARCSIN 

Computer programs are available which generate time histories of the 
rigid body and the mode shapes and frequencies for the body- beam-body 
configuration for "pitch" bending, "roll" bending and "yaw" twisting. 
Since the modes are based on solving explicitly the distributed parameter 
equations (without damping and without kinematic coupling) there is no 
limit to the number of modal characteristic sets that can be generated by 

the program. It will be the analyst’s decision as to how many modes need 
to be considered. 


Laboratory Experiment Description 

The second part of the design challenge is to validate in the 

laboratory, the system performance of the more promising control system 

designs of the first part. The experimental apparatus will consist of a 
dynamic model of the Space Shuttle orbiter with a large antenna reflector 
attached by means of a flexible beam. The dynamic model will be exten- 

sively instrumented and will have attached force and moment generating 
devices for control and for disturbance generation. A single, flexible 
tether will be used to suspend the dynamic model, allowing complete angular 
freedom in yaw, and limited freedom in pitch and roll. An inverted 

position will be used to let the reflector mast to hang so that gravity 
effects on mast bending will be minimized. The dynamics of the laboratory 
model will of necessity be different from the mathematical model discussed 
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earlier. 



Design Challenge, Part One 


For part one of the design challenge, the following mathematical 
problem is addressed. Given the dynamic equations of the Shuttle/antenna 

t t 

configuration, what control policy minimizes the time to slew to a target 
and to stabilize so that the line-of-sight (LOS) error is held, for a time, 
within a specified amount, 5, During the time that the LOS error is 
within 6, the attitude must change 90° to prepare for the next slew 
maneuver. This was previously referred to as the sescondary control task. 
The maximum moment and force generating capability will be limited. Advan- 
tage may be taken of selecting the most suitable initial alignment of the 
Shuttle/antenna about its assigned initial RF axis, line-of-sight. 

Random, broad band-pass disturbances will be applied to the configuration. 
Two proof-mass, force actuators may be positioned anywhere along the beam. 
The design guidelines are summarized below: 

1. The initial line-of-sight error is 20 degrees. 

e(o) = 20 degrees 

2. The initial target direction is straight down. 



3. The initial alignment about the line-of-sight is free to be chosen 
by the designer. Advantage may be taken of the low value of 
moment of inertia in roll. The Shuttle/antenna is at rest 
initially. 

4. The objective is to point the line-of-sight of the antenna and 
stabilize to within 0.02 degree of the target as quickly as 
possible. 

<5 =» 0.02 degree 

3 ?? 
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Cowro! moments can be applied .e ,00 Ha sampling r , t . eo ^ the 
shuttle and reflector bodies of 10,000 ft-lb for each axis. The 
commanded moment for each sals Is limited to 10,000 ft-lb. The 
actual control moment's response to the commanded value Is 
first-order with a time constant of 0.1 second. 

For the rolling moment applied to the shuttle body: 

-10 4 < R, , < io 4 

“* X, 1, command ~ 

M x,l (n + 1} = e " 0,1 \ 2 (n) + (1 - e -0 - 1 ) m (n) 

■ X ’ 1 X, 1 .command' 0 ' 

Equations for other axes and for the reflector body are similar. 
Control forces can be applied at the center of the reflector In 
the X and y directions only. The commanded force In a 
particular direction is limited to 800 lbs. The actual control 

force's response to the commanded value Is first-order with a 
response time of 0.1 second. 

For the side for applied to the reflector body: 


- 80 ° * F Y, command < «» 


F y (n + 1) = e -0 * 1 F (n) + (1 - e" 0 * 1 ) F (n) 

* Y, command'' 

Equations for X-axis are similar. 

Control forces using two proof-mass actuators (each having both 

X and y axes) can be applied at two points on the beam. The 

are limited to r 1 ft, and the masses weight 10 lbs each. 

The actual stroke follows a first-order response to limited 
commanded values* 
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For the X-axis of the proof-mass actuator at S 2 ^ 


-l < A y 0 . < 1 

— X, 2, command — 


A v 0 (n + 1) - e" 0 * 1 A (n) + (1 - e' 0 ' 1 ) A (n) 

X, 2, command 


X,2' 


Equations for other axes and locations are similar. 

8. The inertial attitude direciton cosine matrix for the Shuttle body 
lags in time the actual values by 0.01 second and are made at a 
rate of 100 samples per second. Each element of the direction 

cosine measurement matarix is contaminated by additive, 
uncorrelated Gaussian noise having an rras value of 0.001. The 
noise has zero mean. 



d n (n) 

d 12 (n) 

d 13< n) 

T s .measured^ 11 + ^ T s,true^ n ^ + 

d 21 (n) 

d 22 (n) 

d 23 ( n > 

where: 

d 31 (n) 

d 32 (n) 

d 33 (a) 

E{d ±j (n) } « 0 

E Kj (n)d kL (n) ^ = 0 

for i 

* k or 

j * L 

E Kj (n)d ij (n + k >i - ° 

for k 

* o 


= ( .001 J 2 

for k 

= 0 
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9. The angular velocity measurements for both the Shuttle and 
reflector bodies pass through a first-order filter with 0.05 sec 
time constant and lag in time the actual values by 0.01 second and 
are made at a rate of 100 samples per second. Each rate 
measurement is contaminated by additive, Gaussian, uncorrelated 

noise having an rms value of 0.02 degree per second. The noise 
has zero mean* 

For example: 


l,X>measured (n + ^ ,X, filtered (n) + e i,x (n) 

E K,X ( ") e i,x (n + k) } * 0 for k * 0 


- (.02)2 f or k = 0 

where 


1 ,X, filtered 


— 20 oo + 

1 ,X, filtered 


20 to 


1 ,X, true 


10. Three-axis accelerometers are located on the Shuttle body at the 
base of the mast and on the reflector body at its center. Two- 
axes (X and Y) accelerometers are located at intervals of 
10 feet along the mast. The acceleration measurements pass 
through a first-order filter with a 0.05 second time constant and 
lag in time the actual values by 0.01 second, and are made at a 
rate of 100 samples per second. Each measurement is contaminated 
by Gaussian additive, uncorrelated noise having an rms value of 
0.05 ft/sec^. 
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For example: 


a i ,X, measured^" + ^ a i ,X, filtered^ + T 1,X^ 

E ^ T 1 X (n) T 1 X (n + k) ) = 0 for k * 0 

=* ( .05)^ for k = 0 

where: 

a *i,X, filtered = " 20 3 1 ,X, filtered + 20 “l.X.true 

11. Gaussian, uncorrelated step-like disturbances are applied 
100 times per second to the Shuttle body in the form of 3-axes 
moments, having rms values of 100 ft-lbs. These disturbances 

have zero mean. 

For example: 

E ^ M D,X^ n) M D x < n + = 0 for k * 0 

= ( 100)^ for k = 0 

In summary, the designer's task for part one is to: (1) derive a 

control law for slewing and stabilization, coded in FORTRAN; (2) select an 
initial attitude in preparation for slewing 20 degrees; and (3) select two 
positions for the 2-axes proof-mass actuators. An official system 
performance assessment computer program will be used to establish the time 
required to slew and stabilize the Shuttle/antenna configuration. 
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Design Challenge, Part TVo 


As in part one, the task is to minimize the time to slew and stabilize 
a Shuttle/antenna configuration. The difference is that in part two of the 
design challenge, a physical laboratory model will be used instead of the 
dynamic equations of part one. The constraints on total moment and force 
generation capability will apply to part two, as for part one. Again, the 
analyst may select the initial alignment about the assigned initial RF 
line-of-sight. Disturbances will be injected into the Shuttle/antenna 
model. The designer's task will be similar to that for part one. 

CONCLUDING REMARKS 

A Design Challenge, in two parts, has been offered for the purpose of 
comparing directly different approach to controlling a flexible 
Shuttle/antenna configuration. The first part of the design challenge uses 
only mathematical equations of the vehicle dynamics; the second part uses a 
physical laboratory model of the same configuration. The Spacecraft 
Control Laboratory Experiment (SCOLE) program is being conducted under the 
cognizance of the Spacecraft Control Branch at the NASA Langley Research 
Center. The NASA/ IEEE Design Challenge has the advice and counsel of the 
IEEE-COLSS Subcommittee on Large Space Structures. Workshops will be held 
to enable investigators to compare results of their research. 
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MASS CHARACTERISTICS 











































Figure 2.- Schematic of the effect of bending on the 
1 ine-of-sight pointing error. 
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Roll bending mode ft 4. Frequency= 12.29 Hz 
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4 


4 


* 


Figure 4a.- Plots of normalized roll bending mode shapes 
for SCOLE configuration. 


4 10 



flnPLiTJoc amplitude ftTPHTuot amplitude amplitude amplitude amplitude amplitude 


OF POOR QUALITY 





Pitch bending mode # 1. Frequency = .29 Hi 



Pitch bending mode # 2. Frequency^ 1.65 Hz 



Pitch bending mode # 3. Frequency= 497 Hz 



Pitch bending mode # 4 Frequency= 12.36 Hz 



Pitch bending mode # 5. Frequency= 23.72 Hz 



Pitch bending mode # 6. Frequency= 38.91 Hz 



Pitch bending mode # 7. Frequency^ 57.92 Hz 



Pitch bending mode # 8. Frequency= 80.73 Hz 


Figure 4b.- Plots of normalized pitch bending mode shapes 
for SCOLE configuration. 
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Frequency= 270.69 Hz 


I- 



Torsional mode # 



8. Frequency= 315.80 Hz 


4 


i 


Figure 4c.* Plots of normalized torsional mode shapes for 
SCOLE configuration. 
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The Spacecraft Control Laboratory Experiment Is a facility for the 
Investigation of control techniques for large flexible spacecraft. The control 
problems to be studied are slewing maneuvers and pointing operations. The facility 
Implements the salient characteristics of a flexible satellite with distributed 
sensors and actuators. 

The flexible satellite Is represented by a continuous structure consisting of a 
large mass and Inertia connected to a small mass and Inertia by a slender, flexible 
beam. The structure is suspended by a single cable mounted to a universal joint at 
the system C. G. The sensors for the experiment consist of aircraft quality rate 
sensors and servo-accelerometers. The shuttle attitude will be determined through a 
combination of inertial measurements and optical sensing techniques. Actuators for 
the experiment consist of Control Moment Gyros, reaction wheels, and cold gas 
thrusters. Computational facilities consist of micro-computer-based central 
processing units with appropriate analog interfaces for implementation of the 
primary control system, the attitude estimation algorithm and the CMG steering law. 

Details of the experimental apparatus and the system software are presented in this 
paper* 


ABSTRACT 

A laboratory facility for the study of control laws for large flexible 
spacecraft is described in the following paper. The facility fulfills the 
requirements of the Spacecraft Control Laboratory Experiment (SCOLE) design 
challenge for a laboratory experiment, which will allow slew maneuvers and pointing 
operations. The structural apparatus is described in detail sufficient for 
modelling purposes. The sensor and actuator types and characteristics are described 
so that identification and control algorithms may be designed. The control 
implementation computer and real-time subroutines are also described. 


INTRODUCTION 

A modelling and control design challenge for flexible space structures has been 
presented to the technical community by the NASA and IEEE (ref. 1). The Spacecraft 
Control Laboratory Experiment (SCOLE) was constructed to provide a physical test bed 
or the investigation and validation techniques developed In response to the design 
c a enge. The control problems to be studied are slewing maneuvers and pointing 
operations. The slew Is defined as minimum time maneuver to bring the antenna 
line-of sight (LOS) pointing to within an error limit of the pointing target. The 
second control objective is to rotate about the line of sight and stabilize about 
the new attitude while keeping the LOS error within the bound «. The SCOLE problem 
is defined as two design challenges. The first challenge Is to design control laws, 
using a given set of sensors and actuators, for a mathematical model of a large 
antenna attached to the space shuttle by a long flexible mast. The second challenge 
Is to design and implement the control laws on a structural model of the system in a 
laboratory environment. This report gives preliminary specifications of the 
laboratory apparatus so that interested investigators may begin design and 
simulation for the laboratory experiment. 


a c TJ? lab « ra tory experiment shown in figure 1 attempts to implement the 
definition of the modelling and control design challenge within reasonable limits of 
the 1-g atmospheric environment. The experimental facility exhibits the essential 
SCOLE characteristics of a large mass/inertia (space shuttle model) connected to a 
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small mass/ Inertia (antennae reflector) by a flexible beam. Control sensors and 
actuators are typical of those which the control designer would have to deal with on 
an actual spacecraft. Some trades are made in terms of structure, sensors 
actuators, and computational capability in order to develop the experiment ’ in a 
timely and cost-effective manner. To this end, the basic structure Is made of 
homogeneous, continuous elements. It Is suspended from a steel cable with the 
positive z-axls of the shuttle pointing up, thus minimizing the static bending of 
the antenna mast. The suspension point Is a two-degree-of-freedom gimbal for pitch 
and roll with yaw freedom supplied by the suspension cable. The sensors are 
aircraft quality rate sensors and servo-accelerometers. The shuttle attitude will 
be determined through a combination of inertial measurements and optical sensing 
techniques. 6 


contro1 moments are provided by a pair of two-axis control moment 
gyros (CMG s). Mast-mounted control torques can be applied by a pair of two-axis 
reaction wheels. The reflector-based forces are provided by solenoid-actuated cold- 
air thrusters. Reflector mounted torque devices are a trio of high-authority 
reaction wheels. Computational facilities consist of micro-computer-based central 
processing units with appropriate analog interfaces for implementation of the 
primary control system, the attitude estimation algorithm, and the CMG steering 

law. All of the elements which make up the SCOLE experiment are described in detail 
in the following text. 


The description of the apparatus covers five major groups: The basic 

structural elements are described and pertinent dimensions and structural properties 
are provided. The sensor locations and their dynamic properties are presented. The 
actuator locations and estimated dynamic properties are also given. The mass prop- 
erties of the combined structure, sensor and actuator system are given. Finally 
the computing system and analog interfaces are described. ' 


The contents of this report are considered accurate at the time of 
publication. All of the planned SCOLE components are implemented and are available 
to the user at a raw signal level. However, due to continued refinement of some of 
the components, specific details of the system may change over the life-time of the 
experimental apparatus. 


b 


STRUCTURES 


SCOLE Is comprised of three basic structures, the shuttle, the mast, and 
ector panel. The assembly of these individual components and the global 
reference frame are shown in figure 2. 


The shuttle planform is made from a 13/16-inch steel plate and has overall 

dimensions of 83.8 hv SA.O InrhAC. Tfe fntal * i ~ cm *t « . 

— " v ' a 6" l -*■=> ^vi./ jwuiiQb. me snuttle's 

center-of-mass Is located 3.4 inches below the experiment’s point of suspension, and 
26.8 inches forward of the tail edge (fig. 3). 


The mast is 120 inches long. It is made from stainless steel tubing and weighs 
<*.48 pounds. One-inch thick manifolds are mounted to the mast at each end. The 
assembly of these parts and their dimensions are shown in figure 4. 


The refiector panel is hexagonal in shape, made from welded aluminum tubing 
and weighs 4.76 pounds (fig. 5). It is located 126.6 inches below the SCOLE’s point 
of suspension. The center of the reflector is located at 12.0 inches in the x 
direction and 20.8 inches in the y direction from the end of the mast. 
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The complete system Is suspended from an 11-foot cable attached at the system 
center-of-gravlty via a universal joint. Roll and pitch rotational freedom Is 
provided by pillow-block ball bearings which have an estimated break-out torque of 
0.1 ft-lb. The universal joint is shown in figure 6. It is fixed to the shuttle 
plate, and the system center-of-gravlty is made to coincide with the center-of- 
rotation by means of an adjustable counter balance system. 


SENSORS 

The sensors for the experiment consists of nine servo-accelerometers and two 
3-axis rotational rate sensing units. An optical sensor will provide yaw attitude 
of the shuttle. The power supplies for these sensors are mounted on the shuttle 
plate to minimize the number of large gauge wires which must cross the universal 
joint suspension point. Only a single 115 VAC cable and 33 signal wires cross the 

universal joint. The wires for the sensors are routed on the shuttle and along the 
mas t * 


Accelerometer 8 


All nine accelerometers have 
350 Hz • Linearity is within 0.17 
calibration is presented in figure 
request. 


a frequency response which is nearly flat up 
percent of the full-scale output. A typical 
7. Individual calibrations are available on 


to 


The shuttle-mounted accelerometers shown in figure 8a sense the x y and z 
U "“; /I— «o*or. . 1 . distributed Ire. the su.peSsL^poi* Jo ,ld 
inertial attitude estimation. The locations and sensitive axis are shown in 
tigure 8b. 


at i 1,138 t_ K 0Unted ac «lerometers shown in figure 9a sense x and y acceleration 

at locations about one-third of the mast length from each end. The positions and 
sensing axis of the devices are shown in figure 9b. 

The reflector-mounted accelerometers are shown In figure 10a. They are 
positioned in the center of the reflector below the thrusters and sense the x and y 
accelerations. The coordinates and sensing axis of the devices are shown in 
figure 10b. 


Rate Sensors 


The rotational rate sensors are three-axis, aircraft-quality instruments. The 
r n C A response 18 approximately flat to 1 Hz and -6 db at 10 Hz. Linearity is 
u * ,? € j Ce, J t ^ U H sca l e * A typical calibration is shown in figure 11. The 

rhrKhoId 0%‘VjselT y ‘“ *'* PUCh a * 1S "" d 360 deg/sec - ,or toU - 

three!avl! hU r^?7K U !? ted SenS ° r packa S e « show " ^ «Rure 12a, senses 

! d b y anRular rates of the shuttle plate. Its coordinates and 
sensing axis are presented in figure 12b. 


The mast-mounted rate sensor package, 
angular rates at the reflector end of the 
are presented in figure 13b. 


shown in figure 13a, senses three-axis 
mast. Its coordinates and sensing axis 
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The sensor information required for control system design is summarized in 
Table I. The sensor type is listed in column 2, and its sensed variable is iisted 
in column 3. The analog interface channel is listed next. The coordinates of the 
parts with respect to the universal joint are listed in the next three columns. The 
sensitivities in terms of analog-to-digital converter units are listed next. The 
error list shows the RMS deviation of the rate sensors or the percentage of full- 
scale linearity error for the accelerometers. The linear range of the instruments 
is listed in the next-to-last column. 


OPTICAL SENSOR 

An optical sensor will be provided to determine yaw attitude of the shuttle. 
The optical sensor is a planar photo-diode with appropriate optics mounted on the 
ground. The outputs of the sensor are proportional to the position of an infared 
light source on the shuttle. The sensor data is processed by a dedicated micro- 
controller and then sent to the main CPU over a serial data link for transformation 
to attitude angles. No photographs or calibration data are available for this 
device. 


ACTUATORS 

The actuators consist of both proportional and on-off controllers. Shuttle 
attitude control is provided by a pair of two-axis control moment gyros (CMG’s). 

Mast vibration suppression can be achieved with a pair of orthogonally mounted reac- 
tion wheel actuators positioned at two stations on the mast. Reflector forces are 
provided by four cold gas jets. Reflector torques are provided by three ortho- 
gonally mounted reaction wheels at the end of the mast. As with the sensors, all 
devices are inertial, and the power supplies and amplifiers are mounted on the 
shuttle. Fifteen command signal wires cross the universal joint. All actuators 
were manufactured in house. 


Control Moment Gyros 

The CMG f s each have two gimbals which are equipped with individual direct drive 
DC torque motors. The momentum wheel is mounted in the inner girnbal and driven by 
two permanent magnet DC motors. The nominal operational momentum is about 2.5 ft- 
lb-sec. The girnbal torque motors are driven by current amplifiers so the output 
torque will be proportional to the command voltage sent to the amplifier. The 
girnbal torquers will produce +/- 1.5 ft-lbs at frequencies up to 1kHz. The gimbals 
are Instrumented with tachometers and sine-cosine poteniometers to facilitate 
decoupled control of the shuttle attitude angles. A dedicated computer will be used 
to control the CMG gimbals. Routines will be provided so that users may command 
decoupled shuttle torques or girnbal torque commands. 

The sensitivity calibration curve of a typical girnbal motor is shown in 
figure 14* No other calibration data are available for the CMG's. 

The forward CMG is shown in figure 15a. Note that the outer girnbal is fixed 
and parallel to the pitch axis of the shuttle. The Inner girnbal is nominally orien- 
ted so that the spin axis .of the momentum wheel is parallel to the shuttle z-axis. 
The second CMG is mounted at the rear of the shuttle so that the outer girnbal is 
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^ thC Z T? X f S * ^ i w nCr glmbal 18 nomlnall y oriented so that the rotor 
spin axis is Parallel to the shuttle x-axis. The coordinates and nominal axis of 
actuation of the CMG s are shown in figure 15b. 

Reaction Wheels 

The mast-mounted reaction wheels consist of aluminum disks with inertia of 

about 0.00027 lb-ft-se c 2 mounted directly on the drive shaft of a 20 oz-in 
permanent magnet DC motor. The motors are powered by high bandwidth current 
A torque sensitivity plot is presented in figure 16. No other 
d TK 3 are avallable * A typical reaction wheel assembly is shown in 

"gu” !Jb! '"° * CtU “ 0,: 1 °" tl0nS a " d thelr aala « f «. 

ss a siip ™ ^ 

caM h?m!°^ r !K are powered by high bandwidth current amplifiers. The torque 

oJCi «?*h ef 5 ! C6S 18 estlmated t0 * ab out 50 oz-in. No sensity^lot or 
J al j bra 5 lon data la Presently available. The three-axis reaction wheel 

Sis. ;„*5n; r e 18 “r te i8a - Tha actuator i °' a “- >aa a " j ^ <* 


Thrusters 

lets Th Thr?^°\ f ° rCeS ° n the J eflector are Provided by solenoid actuated cold gas 
jets. The thrusters are mounted in the center of the reflector and act in the x-v 

n^r* / a J f 8 arG f UppllGd by a compressed air tank mounted on the shuttle. The 

l™ Tllt'tr C r elS the maSt t0 the solenoid -nifold, which gat^ tie 

* lo J b ? t "** n the regulated supply tank and the thrusters as shown in figure 19. 

Vll tr isln J tiatad by opening the solenoid with a discrete command. The rise time 

tlon of t^ hr SC r i a f 1<>n °t 18 ^ ln fl * ure 20 * The -Snitude and d«H 

of the thrust before the air supply is depleted at 60-psi nozzle pressure is 

Table II. ^ ® perCinent data from figures 20 and 21 are tabulated ln 


are sh^^n^reMb^ ^ Their l0Catl °" a " d a * is Nation 

lnformatlon required for control system design is summarized in 
Ina L £ J 8C h tUat °? tyPe and dlrectlon of action* are listed in Column 2 ft. 

r U J t channels are listed in Column 3. The coordinates of the devices 

of d!v(j ? i he , ne three columns - The sensitivities of the actuators in terms 
of digital-to-analog converter units are shown in Column 7. An estimate of the 
tbr„st.r RHS deviation. exhibited It figure ,9 1. bres.nted „ e ro d a! So other 

ZllUtl:rjZi UiU - Iha ™ aXl "'” af tha a >- a ‘ a » a ““ a ^- a * in Che 
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MASS PROPERTIES 


4 


K 


form the ^0LE i apparlt^J R ire°LtaLg^d 1 in S Table e iV Of ® < * ulpment * whlch collectively 
the point of suspension * to* the Approximate I” * eaSUred fr °" 

coordinate^are ^ h° f ^ ^ ^ ^ nt ‘ ^ 

The remaining columns are tte mass moments t^als" f* 

the complete system are presented on the bottom row. lnertia * The totals for 


oldiui 


- irzzziiz ?^r:rrr baMd 

^ z'zz rr u ^ z ™ r ss 

has 12 serial ports and 1 parallel port Te^mf C ® ti<>n8 P r °8 ramm ing* The computer 

b n &“,*“?■ - 8 

CMC control software r.qul^ f« Jta 2 «S»r. 23. The 
controls dsslqn.r who will b. oper.tln* .„ t* Cpd/co^,,"'^ t0 tha 

.-,u« r ”t“^i'“*^“; 1 ‘2 d ? n SSSiiT'ST ‘ 1 settln ? th * "*“•» 

listed below. PP di A * Th inost common ly used routines are 


For accessing the analog devices: 

getadc - read the analog-to-dlgital converters 
set ac set the digital— to-analog converters 
thrust - set the cold-gas thrusters. 

To control the sampling interval: 

rtime - sets the sample period marks the beginning 
of a real-time loop. 

£ k :- 

sr^rsu-..-,: 

tJ.CF'-% a ur:^: o 'L„ T ^^r r :“*;c:.r;““L 8 :“ u np“ i1 of 

commands are complete, the routine must return to the top of the real-time 
loop and once again call rtime. If this occurs before the end of the 

lTll U f i n ^ rVa \ th ? “ me "° ut condition will be inhibited and rtime will 

for the next rising edge of the sample Interval clock and then return 
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to the calling program. If the user computations take longer than the 
sample interval, a time-out condition will be signaled to the operator 
when rtime is called. The user may choose to Ignore the condition and 
continue or may take specific steps to alleviate the condition. 

The procedure for logging on to the computer is as follows: 

Set communication parameters to 1200 baud, 7 bit, even parity. 

Dial in to the Langley data communication switching system at 
804-865-4037. When connected, type a carriage return. 

To the system prompt ’‘ENTER RESOURCE CODE*' type “acrl." 

Wait until “GO" and the “name:** prompt appear on the screen. 

Type in your log-in Information. All investigators will be given a three- 
letter log-in name (usually the university affiliation). 

Some useful system commands are listed below: 

To transfer a file to the experiment computer from a smart terminal, type 
cat> filename <cr> and then enable the upload function of the local 
terminal. When the upload is complete, type <cntl>d. 

To transfer a file from the experiment computer to a smart terminal, type 
cat filename then enable the download function of the local terminal and 
type a <cr>. 

To list the contents of a directory, type 1 <cr>. 

To look at a file, type p filename. 

To compile a FORTRAN program and link with real-time system commands and 
TCS graphics, type frt filename. 

Note: filename must have the extension .for. The executable code will be 

under filename without the .for extension. To run, simply type filename 
without any extension. 

To list the system commands, type 1 /bin. to get a description of any 
command, type describe cownd. 


System user guides will be available upon request from the Spacecraft Control 
Branch, M/S 161, NASA Langley Research Center, Hampton, VA, 23665-5225. Other 
details for operating in a real-time mode will be provided at the time of 
implementation. 


CONCLUDING REMARKS 

The SCOLE laboratory facility is an experimental apparatus which permits 
ground-based investigation of identification and control algorithms for large space 
structures. The facility exhibits structural dynamics similar to those expected on 
the large satellites. The sensors and actuators are typical of those, which may be 
used on an operational satellite. The computational system Is reasonably sized with 
current technology processors and permits ready access to the facility for 
interested investigators. 
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The description of the structural assembly, the sensor and actuator 
configuration, and software provided in this paper should be sufficient for 
Investigators to begin designing identification and control algorithms for 


SCOLE 
the SCOLE 
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APPENDIX 


REAL-TIME SYSTEM SUBROUTINES 

Analog I/O system command. 


NAME: 

getadc Samples the analog— to— digital converters. 


COMMAND: 

getadc [ — s ] 


DESCRIPTION: 


c£ann^r a r d * S ^ed to sample the analog- to-digltal converters and display selected 
channels at the terminal. The +/- 10.0 volt Input range is scaled to +/- 1 n tl,ll 

with Jh ngl€ blt . ls wor ^ h •° 0049 unit8 ‘ The channels to be displayed are selected 
orint n " 8 ° ptl0 "‘ , Tbls °P tlon displays a menu which allows the user to set the 
print flags for the individual ADC channels. Specific choices are: 


1) turn on all print flags, 

2) turn off all print flags, 

3) turn on a range of print flags, 

4) turn off a range of print flags, 

5) display current print flags, 

6) save current print flags. 


If the command Is executed 
used to selectively display 


without the -s option, the last set of 
the ADC channels. 


print 


flags 


is 


USES: 

getadc. flags 


DIRECTORY: 

/bin 


SOURCE: 

/ usr/ csc/ele/getadc.c 
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Analog 1/0 system command. 


NAME: 

getadc Samples the analog-to~dlgital converters. 


COMMANDS 

getadc [-s] 


DESCRIPTION: 


I5in«!?! ma r d .i 8 “ 8ed /° ® a “Ple the analog-to-digital converters and display selected 
i i* K 4 * * erminal * The + /~ 10 *° volt input range is scaled to +/- 1.0 units 
? * H ng blt , ls worth *00049 units. The channels to be displayed are selected 
with the -s option. This option displays a menu which allows the user to set the 
print flags for the individual ADC channels. Specific choices are: 


1) turn on all print flags, 

2) turn off all print flags, 

3) turn on a range of print flags, 

4) turn off a range of print flags, 

5) display current print flags, 

6) save current print flags. 


If the command is executed without the -s option, 
used to selectively display the ADC channels. 


the last set of print flags is 


USES: 

getadc. flags 


DIRECTORY: 

/bln 


SOURCE: 

/usr/csc/ele/getadc.c 
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Analog I/O system command. 


NAME: 

setdac Sets the digital— to-analog converters. 


COMMAND: 


setdac [— s] [-0] 


DESCRIPTION: 


Ihann^l! Jk 1 */ Se ? n the Voltage on a ran «« of digital-to-analog output 

channels. The +/- 1 .0 unit output range is scaled to +/- 10.0 volts so a single 

unit is worth .0049 volts. The channels to be set are selected by executing the 
command with the -s option. This option displays a request for the range of 
channeis to be set, and then queries for individual channel values in terms of 

units. All DAC channels may be set to zero by executing the command with the ••-0• 
option. 


USES: 

Nothing. 


DIRECTORY: 

/bin 


SOURCE: 

/usr/ csc/ele/setdac.c 
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UNOS system command. 


NAME: 

* 


Iter. Terminal emulator. 


1 


COMMAND: 

term [-si [-S] 


DESCRIPTION: 


This command connects the user terminal to the Langley central data communication 
switch at 1200 baud. This is a dumb terminal emulator which provides rudimentary 
file transfer capabilities. No attempt is made to emulate control codes of any 
particular terminal for editing purposes. 

The emulator commands are as follow: 

<? Return to UNOS (operating system.) 

! To download a file. 

To upload a file. 

To excute a system command. 

? For help. 

The options are: 

•s 300 baud 

-S 900 baud 

Upload means to transfer a file from the Charles River computer to the remote 
computer. The remote computer must have some mechanism for receiving the text. 

Download means to transfer a file from the remote computer to the Charles River 
computer. No attempt is made to check for existence of the receiving file name 
before saving the downloaded file. 


USES: 

/ doc/ cmds/lterm.help 


DIRECTORY : 

/bin 


r 


SOURCE: 

/jpw/ltenn.c 


4 
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Analog I/O system command. 


NAME: 

8751_test Test the serial communication link to the 8751 boards. 


COMMAND: 


8751 test 


DESCRIPTION: 

C r and verlflcatl0n and callbra tl°" of the 8751 micro-controller 

ce over the RS-232 serial ports. The command queries for voltages to be 
output by the digital-to-analog converters on the 8751 boards. The data input is in 
terms of units with 2047 equal to 9.9951 volts and -2048 equal to -10.0000 volts. 

ILirhi* 1 * , ? S ° nly ° ne V3lue ’ 311 actlve boards are sent that value. 

Otherwise, individual values are sent. 


USES: 

raotint( ) 
motsub( ) 


DIRECTORY: 

/bin 


SOURCE: 

/usr/rdb/8751 test. for 


i 


* 
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Analog I/O system subroutine 


NAME: 

getadcQ ( C callable ) Sample a range of analog-to-digital converters. 


CALL: 

int error, flrst_adc, last_adc; 
float adc_data_pointer; 

Int getadc( first_adc, last__adc, &adc_data_pointer) 
error - getadc(f irst_adc, last_adc, &adcjdata_pointer) 


DESCRIPTION: 

This subroutine samples a range of analog-to-digital converters. The +/- 10.0 volt 
Input range Is scaled to +/- 1.0 units so a single bit is worth .00049 units. The 
arguments are: 

first_adc (int) First converter to be sampled (numbering starts from 

zero.) 

l ast _ ac * c (Int) Last converter to be sampled (maximum is 63.) 

&adc_data_pointer (*) Starting location for storing sample data. Data are 

floating point values with a range of +/- 1.0. 


RETURNS: 

error * 0 indicates valid transfer, 
error * indicates bad range. 


USES: 

Nothing 


LIBRARY: None 

SOURCE: 

/usr/csc/ele/getadc_c.c 
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Analog I/O system subroutine. 


NAME: 

setdacO 


( C callable ) Set a range of dlgital-to-analog con vet era. 


CALL: 

Int error, first dac, last dac; 
float dac data pointer; 

lnt setdac( first dac, last dac, & dac data pointer ) 


DESCRIPTION: 


This subroutine sets s range of dlgit.l-to-anelog converters. The +/- ! n unit 
argument^are; 8 “ +/ ‘ ,0 -° ''° US “ 8 «“ * .O^'CSlt”!' The 


first dac (lnt) 


First converter to be set (numbering starts from 
zero.) 


last dac 


( int) 


Last converter to be set (maximum Is 7.) 


&dac_data_pointer (*) 


Starting location DAC data. Data are floating point 
values with a range of +/- 1.0. 


RETURNS: 

error = 0 Indicates valid transfer, 
error =-l Indicates bad range. 

error > 0 Indicates "error” number of data words out of range 

USES: 

Nothing, 

LIBRARY: 

None. 


SOURCE: 

/ usr/ csc/ele/ set dac c.c 
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Analog I/O system subroutine. 


NAME: 

getadcO ( Fortran callable ) Sample a range of analog-to-digital 

converters. 

CALL: 

integer error, getadc 

error - getadc( f irst_adc, iast_adc, adc_data_array) 

DESCRIPTION: 

This subroutine samples a range of analog-to-digital converters. The +/- 10.0 volt 
Input range is scaled to +/- 1.0 units so a single bit is worth .00049 units. The 
arguments are: 

adc (integer) First converter to be sampled (numbering starts from 

zero.) 

last_adc (integer) Last converter to be sampled (maximum is 63.) 

adc_data_array (real) Starting location for storing sample data. Data are 

floating point values with a range of +/- l.o. 

RETURNS : 

error = 0 indicates valid transfer, 
error *-l indicates bad range. 

USES: 

getadc^w. j . 

LIRRARY: 

/iib/ acrl__rt_JLib_f . j 

SOURCE: 

/ usr/csc/ele/getadc_f .c 
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Analog I/O system subroutine. 


setdacO ( Fortran callable ) Set a range of digital— to— analog 

converters. 


CALL: 

integer error, setdac 

error - setdac( flrst_dac, last_dac, dac_data arrray ) 


DESCRIPTION: 


This subroutine sets a range 
output range is scaled to +/- 
arguments are: 


of digltal-to-analog converters. The +/- 1.0 unit 
10.0 volts so a single unit is worth .0049 volts. 


The 


first_dac (integer) First converter to be set (numbering starts 

from zero.) 

last_dac (Integer) Last converter to be set (maximum is 7.) 

dac_data_array (real) Starting location DAC data. Data are floating point 

values with a range of +/- 1.0. 


RETURNS: 

error = 0 Indicates valid transfer, 
error *-l Indicates bad range. 

error > 0 Indicates error number of data words out of range 

USES: 

setdac__w, j 


LIBRARY: 

/lib/acrl_rt_J.ib_f . j 


SOURCE: 

/usr/csc/ele/ setdac f.c 
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Analog I/O system subroutine 


NAME: 

votsubO ( Fortran callable ) Send scaled voltages to the 8751 

micro-controllers which in turn set 
individual DACs. 


CALL: 

integer torque 

call motsub ( torque ) 


DESCRIPTION: 

This subroutine sends the motor torque command data to the 8751 micro-controller 
boards which in turn load the data into the digital— to— analog converters. 

torque(6) (integer) Array of dimension 6 which contains the seated data 

to be output on the 8751 DAC’s. The range of the data 
is +2047 for +9.9951 volt output to -2048 for -10.0000 
volt output. 


REQUIRES: 

Call to motint. 


RETURNS: 

Nothing . 


USES: 

motsub w.j 


LIBRARY: 

/ lib/acrl__rt_lib_f . j 


SOURCE: 

/jpw/8751 COM/torsub.c 
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Analog I/O system subroutine. 

NAME: 

finlshO ( Fortran callable ) Close serial ports to 8751s. 

CALL: 

finishO 

DESCRIPTION: 

This subroutine closes the serial communication lines to the 8751s. There are no 
arguments. 

RETURNS : 

Nothing. 

USES: 

finish__w.j . 

LIBRARY: 

/lib/acrl_rt_lib_f . j 

SOURCE: 

/jpv/8751COM/torsub.c 
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Analog I/O system subroutines. 


NAME: 

thrust ( Fortran callable ) Set the discrete ports to activate the 

thrusters. 


CALL: 


integer*2 thrust * thrust ( x,y ) 


DESCRIPTION: 

This subroutine sets the states of the four discrete outputs on the Parallel 
Interface/Timer. It was designed for the thrusters on the SCOLE facility. The 
arguments can have one of three values: 1, 0, or -1 corresponding to positive, 

none, and negative thrust respectively. 

The arguments are: 

x (integer) State of x thruster. 

y (Integer) State of y thruster. 


REQUIRES: 

Nothing. 


RETURNS : 

Nothing . 


USES: 

thrust_w. j 


LIBRARY: 

/lib/acrl rt lib f . 


SOURCE: 

/jpw/TIMER/ pitasc. j 
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Analog I/O system subroutine. 


NAME: 

rtl«e() (fortran callable) Mark the start of the real-time loop and 

change the sample time Interval If required. 

CALL: 

call rtime( tau, flag ,k) 

DESCRIPTION: 


mapped tlmer ls used to control the timing of real-time 

gene rate "a start ^ arles , River Computer. The programmable ^lock is required to 

process The S ? V a at ° p pulse for each sampling interval of the control 

process. The maximum interval Is eighty-five seconds and the minimum interval Is 5 

micro-seconds. The timer is a Motorola M68230 Parallel Interfac Timer (PI/T) uhlrh 

provides versltlle double buffered parallel Interfaces and 74— hit nm ki 

for M68000 systems. interfaces and 24-bit programmable timer 

arguments 'a re^ ^ ^ Sh ° Uld * ”*** JuSt lnslde the real time loop. The 


tau 

flag 

k 


(real) is the sample period, 

(loglcai) is the indicator to either maintain the same value of tau or 
pass in a new value, 

(integer) is the timout parameter. If k is returned from rtime 

a°aiouJ n has i* thl8 lndl * ate * a nor ®al return. If k ls returned a 0, 
a timout has occurred, and appropriate action should be taken. The 
user must supply his/her own timout procedure. 


REQUIRES: 

Nothing. 


RETURNS: 

k 


LIBRARY: 

/lib/acrl_rt_lib_f . j 
Add " /j pw/rtirae.obj 


to FORTRAN compile command. 


USES: 

traset ( ) 
init() 
rtwate( ) 
cktim( ) 
pint_w. j 

SOURCE: 

/j pw/TIMER/pint.c 
/ jpw/TIMER/pint__w.m 
/jpw/ rtime. for 
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SENSOR PARAMETERS FOR CONTROL SYSTEM DESIGN 



























































































ZABLE II. TYPICAL THRUSTER DATA 
TABLE II. TYPICAL THRUSTER DATA FDR SIXTY PSI NOZZLE PRESSURE. 


Peak Thrust 
Steady State Thrust 
Rise Time 
Thrust Duration 


0.641 lb 
0.32 lb 
0.032 seconds 
24 seconds 


#> 
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TABLE HI. ACTUATOR PARAMETERS FOR CONTROL SYSTEM DESIGN 














































































TABLE IV. SCOLE APPARATUS MASS PROPERTIES AND COMPONENT LOCATIONS 



























































































SCOLE APPARATUS HASS PROPERTIES AND COMPONENT 



4 - 4-0 











































































































































































































TABLE IV. SCOLE APPARATUS MASS PROPERTIES AND COMPONENT LOCATIONS (CONT'D) 



0.0 - 3.4 501.73 2458.5 0.0 - 1643.3 1402.8 6467.7 8510.2 0.0 254 







































































































































TABLE IV. SCOLE APPARATD 
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The SCOLE experiment 


Figure 1 


apparatus 
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1.68 lb 


Figure 4, Mast and manifold 


assembly 



Bored cylinder 
4.00 outer diameter 
t = 0.25 depth = 1.50 
Top surface t = 0. 25 — \ 


Tubing A1 
t = 0.0625 


O.750 outer diameter 



Figure 5. Ref le 
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Figure 6* Universal joint suspension point* * 
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Figure 7. Typical accelerometer calibration 



Figure 8a. Shuttle-mounted accelerometers. *huttle-raounted accelerometers 



Figure 9a. Mast-mounted accelerometers. 



9b - M sensing .xl. 

Of mast-mounted accelerometers. 
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Figure 10b. Coordinates and sensing axis of 

Figure 10a. Reflector— mounted accelerometers. 

reflector mounted accelerometers. 
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Figure 11. Typical rate sensor calibration. 
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Figure 12a. Shuttle-mounted, three-axis 
rate sensor. 



Figure , 2b . Coordinates and e.„. lng „„ of 
shuttle-mounted rate sensor. 



Figure I3a. Mast-end-mounted, three-axis 
rate sensor. 
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Figure 13b. Coordinate, end „ ls of 

mast-end-mounted rate sensor* 
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Figure 14. Typical CM G giuibal 
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Figure 16. Typical reaction wheel torque sensitivity curve. 



Figure 17a. Mast-mounted reaction wheels. 


Figure 17b. Coordinates of mast-mount f ed 
reaction wheels. 
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Figure 18b. Coordinate and axis of actuation 

for mast-end mounted reaction 
wheels. 


Figure 18a. Mast-end mounted reaction wheels. 



Figure 20. Thruster startup and transients 



Figure 19. Thruster air supply schematic. 



Figure 21. Thrust magnitude and duration. 
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Figure 22a. Reflector-mounted thrusters. 



Figure 22b. 


Coord!..!., „ f refl«ctor-mount«d thn.st.rs 
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Figure 23. SCOLE computer interfaces 
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"AWARDS" 


A. V. Balakrishnan was awarded full membership into the 
ORDER of the INFINITE ORDER SYSTEMS. 

Bob Skelton was awarded full membership into the ORDER 
of the FINITE ORDER SYSTEMS. 

R. K. Miller was presented and award for SMOOTHING by 
BANG-BANG Control. 

Peter Bainum was presented an award for studying the 
SIGNIFICANCE of nearly INSIGNIFICANT effects. 

Gene Lin was presented an award for the SLEWING of 
MONSTROUS structures. 

Suresh joshi was presented an award for LOOKING at the 
LINE-of-SIGHT. 


Howard Kaufman was presented an award for being a 
MODAL MODEL FOLLOWER. 

Mike Fi sher was membership into the order of 

the MYSTERIOUS, for his studies into the MYSTERIES of 
actuator placement. 
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Summary of Panel Discussions 

1. The mathematical problem statement is incomplete in 
the sense that "rigidizaUon" of the equations of motion 
does not yield the rigid body equations. Larry Taylor and 
Yogi Kakad intend to derive the missing kinematic terms 
in the partial differential equations of the flexible beam. 

2. An accurate "proof model" of the SCOLE is needed for 
comparing approximate model responses and to evaluate 
the various control laws. Larry Taylor will see that such 
a model is made available as soon as a suitable model 
exists. The modeling difficulties have been more 
troublesome than expected. 

3. The use of the term "modes" has been used loosely to 
include the admissible functions which are often tnndes 
under condition of no damping and no control. Also, the 
full model is nonlinear and consequently does not exhibit 
modes but mode-like characteristics. 

4. It will be necessary to employ relatively low-order 
models for control synthesis, state estimation and on-line 
control. One must use caution, however, lest the 
designer forgets that the full model or the actual 

experimental apparatus is nonlinear and has infinite 
order. 

5- Caution is caned for in employing Kalman filter 
techniques for time-varying, nonlinear dynamics of 
infinite order, such as the SCOLE problem. 
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